SCHOOL 

OF 

OPERATIONS  RESEARCH 
AND 

INDUSTRIAL  ENGINEERING 


docum#nt 


COLLEGE  OF  ENGINEERING 
CORNELL  UNIVERSITY 

ITHACA,  NEW  YORK  14853 


g?,p>90V<W  "woo  31IJ  3Qfl 


SCHOOL  OF  OPERATIONS  RESEARCH 
AND  INDUSTRIAL  ENGINEERING 
COLLEGE  OF  ENGINEERING 
CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK 


TECHNICAL  REPORT  NO.  394 

December  1978 


AN  ANALYSIS  OF  ONE  WAREHOUSE, 

N RETAILER  PRODUCTION  INVENTORY  SYSTEMS 

by 

Howard  Mark  Singer 


D D C 

ir?rR>r?nrmr?|m 

FEB  18  1SI9 

lUTTE 


riiia  document  hau  been  apT>r^vpd 
for  public  release  and  ealo;  its 
distribution  is  unlimited. 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under 
contract  N00014-75-C-1172 , Task  NR  042-335. 


79  02  06  006 


In  this  thesis  we  examine  the  one  warehouse,  N retailer  production 
inventory  system.  We  present  properties  of  optimal  operating  policies 
and  methods  to  determine  various  operating  policies  given  a fixed  cost 
for  set-up  and  an  inventory  carrying  charge  at  each  facility.  We 
assume  that  the  external  demands  on  this  system  occur  either  at  a known 
continuous  rate  that  is  stationary  over  an  infinite  time  horizon,  or 
at  a known  rate  that  may  vary  in  each  of  a finite  number  of  periods . 

We  will  refer  to  the  former  case  as  a continuous  one  warehouse , N 
retailer  problem  and  the  latter  as  a dynamic  demand  problem. 

A one  warehouse,  N retailer  system  is  a special  case  of  the  more 
general  arborescent  production  inventory  system.  We  begin  by  examining 
the  previous  research  on  systems  with  this  arborescent  structure  as 
well  as  reviewing  the  literature  dealing  with  the  serial  and  assembly 
multi-echelon  production  inventory  systems. 

The  continuous  one  warehouse , N retailer  problem  is  then  examined 
in  detail.  The  basic  model  is  introduced  as  well  as  some  previously  solved 
special  cases.  Several  basic  production  policies  that  have  been  suggested 
for  this  system  are  reviewed  and  properties  of  an  overall  optimal  solution 
are  discussed.  We  develop  optimal  and  heuristic  algorithms  to  determine 
the  values  of  the  parameters  in  "single  cycle"  policies  (Schwarz,  1973) 


and  make  comparisons  based  on  both  the  quality  of  solutions  obtained 
and  the  computational  effort.  It  has  been  conjectured  by  Graves  and 
Schwarz  (1977)  that  these  single  cycle  policies  are  optimal  for  certain 
larger  classes  of  production  plans  and  we  demonstrate  that  this  is  not 
the  case.  Finally,  we  aiscuss  the  class  of  multiple  cycling  policies 
and  their  relation  to  optimality. 

We  then  explore  the  dynamic  one  warehouse,  N retailer  problem. 

The  basic  model  is  introduced  and  both  the  Wagner-Whitin  (1958)  and 
the  Zangwill  (1966)  algorithms  are  reviewed.  Veinott  (1969)  has 
suggested  the  extension  of  Zangwill 's  algorithm  to  the  case  of  N 
retailers,  and  we  present  this  extension  along  with  computational  simpli- 
fications. Easily  computed  upper  and  lower  bounds  on  the  optimal  solution 
value  are  developed.  Exact  solution  techniques  for  several  special  cases 
for  the  cost  structure  are  also  discussed.  Finally,  areas  for  future 
investigation  are  suggested. 
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CHAPTER  1 


INTRODUCTION 

1.1  General  Introduction 

In  this  thesis  we  examine  the  one  warehouse,  N retailer  production 
inventory  system.  We  present  properties  of  optimal  operating  policies 
and  methods  to  determine  various  operating  policies  given  a fixed  cost 
for  set-up  and  an  inventory  carrying  charge  at  each  facility.  We  assume 
that  the  external  demands  on  this  system  occur  either  at  a known  con- 
tinuous rate  that  is  stationary  over  an  infinite  time  horizon,  or  at  a 
known  rate  that  may  vary  in  each  of  a finite  number  of  periods . We 
will  refer  to  the  former  case  as  a continuous  one  warehouse , N retailer 
problem  and  the  latter  as  a dynamic  demand  problem. 

A one  warehouse,  N retailer  system  is  a special  case  of  the  more 
general  arborescent  production  inventory  system.  In  the  next  section  we 
examine  the  previous  research  on  systems  with  this  arborescent  structure 
as  well  as  review  the  literature  dealing  with  the  serial  and  assembly 
multi-echelon  production  inventory  systems.  Particular  emphasis  will  be 
placed  on  the  assumptions  in  these  works  that  are  common  to  the  research 
in  this  thesis. 

In  Chapter  2 the  continuous  one  warehouse , N retailer  problem  is 
examined  in  detail.  The  basic  model  is  introduced  as  well  as  some 
previously  solved  special  cases.  Several  basic  production  policies  that 
have  been  suggested  for  this  system  are  reviewed  and  properties  of  an 
overall  optimal  solution  are  discussed.  We  develop  optimal  and  heuristic 
algorithms  to  determine  the  values  of  the  parameters  in  "single  cycle" 
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policies  (Schwarz,  1973)  and  make  comparisons  based  on  both  the  quality 
of  the  solutions  obtained  and  the  computational  effort.  It  has  been 
conjectured  by  Graves  and  Schwarz  (1977)  that  these  single  cycle  policies 
are  optimal  for  certain  larger  classes  of  production  plans  and  we  demon- 
strate that  this  is  not  the  case.  Finally,  we  discuss  the  class  of 
multiple  cycling  policies  and  their  relation  to  optimality. 

In  Chapter  3 we  explore  the  dynamic  one  warehouse , N retailer 
problem.  The  basic  model  is  introduced  and  both  the  Wagner-Whitin 
(1958)  and  the  Zangwill  (1969)  algorithms  are  reviewed.  Veinott  (1969) 
has  suggested  the  extension  of  Zangwill' s algorithm  to  the  case  of  N 
retailers,  and  we  present  this  extension  along  with  computational  simpli- 
fications. Easily  computed  upper  and  lower  bounds  on  the  optimal  solution 
value  are  developed.  Finally,  exact  solution  techniques  for  several 
special  cases  for  the  cost  structure  are  discussed. 

In  Chapter  4 we  summarize  the  major  results.  In  addition,  we  suggest 
areas  for  future  investigation. 

1.2  Serial,  Assembly  and  Arborescent  Production  Inventory  Systems 

A multi-echelon  production  inventory  system  is  one  in  which  products 
are  partially  processed  at  certain  facilities.  Each  subsequent  facility 
performs  additional  processing  until  the  product,  or  products,  is  available 
to  satisfy  external  demands.  Work  in  process  inventories  are  permitted 
between  the  facilities.  The  problem  is  to  determine  an  operating  policy 
that  meets  the  external  demands  while  minimizing  production  and  inventory 
costs . 
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All  the  research  on  multi-echelon  production  inventory  systems  is 
based,  at  least  in  part,  on  the  methods  used  in  explicitly  solving  the 
single  facility  problem.  Single  facility  problems  have  been  studied 
extensively  and  excellent  summaries  of  the  results  on  these  systems 
can  be  found  in  the  books  by  Hadley  and  Whitin  (1963),  and  Johnson  and 
Montgomery  (1974).  For  the  case  of  dynamic  demands,  the  classic  paper 
is  by  Wagner  and  Whitin  (1958).  Their  results  have  been  generalized  many 
times  to  include  multiple  set-up  costs  (Lippman,  1969),  piecewise  concave 
costs  (Love,  1973)  and  backlogging  (Zangwill,  1966).  One  generalization 
of  the  basic  single  facility  problem  that  has  been  extensively  researched 
is  the  capacity  constrained  problem  (Florian  and  Klein,  1971;  Lambrecht 
and  Vander  Eecken,  1978;  Baker,  Dixon,  Magazine  and  Silver,  1978).  The 
Wagner-Whitin  single  facility  algorithm  and  a capacity  constrained  single 
facility  algorithm  can  be  used  to  determine  optimal  policies  in  one 
warehouse,  N retailer  systems  under  certain  cost  restrictions.  (These 
results  will  be  developed  in  Chapter  3.) 

One  can  envision  many  configurations  for  the  facilities  in  a multi- 
echelon system;  but,  three  particular  structures  have  been  extensively 
studied  in  the  past.  These  three  structures  are  the  serial,  assembly 
and  arborescent  systems  (Figure  1.1).  We  will  discuss  each  of  these 
structures  individually  and  summarize  the  relevant  research. 

In  a serial  multi-echelon  system,  each  facility  has  only  one  immediate 
predecessor  and  one  immediate  successor.  Work  in  process  inventory  flows 
from  one  facility  to  the  next  and  at  each  facility  additional  processing 
is  performed.  External  demands  occur  only  at  the  final  facility.  A 
single  conveyor  belt  assembly  line  with  work  stations  along  the  route 
of  the  belt  is  an  example  of  a serial  multi-echelon  system. 
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(a)  SERIAL: 


EXTERNAL  DEMAND 


(b)  ASSEMBLY: 


(c)  ARBORESCENT: 


EXTERNAL  DEMANDS 


FIGURE  1.1  MULTI-ECHELON  SYSTEM  STRUCTURES 
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The  research  on  serial  production  inventory  systems  with  continuous 
deterministic  demands  includes  two  restrictive  assumptions.  Operating 
policies  are  chosen  from  among  the  policies  in  which  there  is  a constant 
lot  size  at  each  stage  and  the  lot  size  at  a stage  is  an  integer  multiple 
of  the  lot  size  at  the  succeeding  stage.  Although  this  type  of  policy 
may  be  optimal  for  concave  cost  cases , there  is  no  reason  to  assume  it 
is  optimal  in  a multistage  system  having  constraints  imposed  by  fixed 
production  rates  and  lead  times. 

In  addition  to  the  assumptions  discussed  above,  Taha  and  Skeith  (1970) 
also  assur  s that  there  is  a delay  between  production  and  use  of  a lot  and 
that  backorders  are  permitted  at  the  final  stage.  No  explicit  solution 
technique  is  given  to  determine  the  integers  multiplying  the  lot  sizes 
that  specify  the  policy.  Jensen  and  Khan  (1972)  examine  a similar  con- 
figuration but  without  the  assumption  of  a delay  between  production  and 
use  of  a lot.  A discrete  dynamic  programming  formulation  to  determine 
the  optimal  cycle  length  is  developed  for  what  is  essentially  a continuous 
parameter  problem.  Johnson  (1972)  provides  expressions  for  the  average 
in-process  inventory  levels  with  or  without  delays  following  production. 

A direct  solution  is  obtained  for  the  two  stage  system  and  an  iterative 
procedure  is  developed  for  the  three  stage  system. 

For  the  case  of  dynamic  demand,  the  classic  reference  is  by  Clark 
and  Scarf  (1960).  They  introduce  the  concept  of  "echelon  stock"  for 
computational  simplicity  and  prove  the  optimality  of  a modified  (s,S) 
policy  for  a serial  system  in  which  the  demand  in  each  period  is  determined 
by  an  observation  from  a known  distribution.  The  echelon  stock  at 
facility  i is  inventory  that  has  been  processed  by  facility  i and  is 
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awaiting  processing  at  any  succeeding  facility.  That  is,  the  facility  i 
echelon  stock  is  inventory  that  has  passed  through  facility  i but  has 
not  yet  been  used  to  satisfy  external  demands.  An  (s,S)  policy  is 
being  followed  if,  at  review,  the  stock  on  hand  plus  on  order,  i,  is 
less  than  s,  then  S-i  units  are  ordered;  otherwise,  no  ordering  is 
done.  A more  detailed  discussion  of  (s,S)  inventory  models  can  be 
found  in  the  book  by  Tijms  (1976). 

For  the  serial  system  with  deterministic  dynamic  demand,  the  basic 
references  are  by  Zangwill  (1966,1968,1969).  He  extends  the  single 
facility  Wagner-Whitin  model  (1958)  to  the  serial  system  with  backlogging 
permitted  at  the  final  stage.  Of  particular  interest  are  the  representation 
of  the  problem  as  a minimum  cost  flow  in  a single  source,  concave  cost 
network  and  the  characterization  of  the  solution  by  exact  requirements , 
that  is , the  production  quantities  are  always  the  sum  of  the  requirements 
for  a certain  number  of  future  consecutive  periods.  The  network  repre- 
sentation and  exact  requirements  will  be  discussed  in  more  detail  in 
Chapter  3 . 

As  in  the  single  facility  case,  there  have  been  many  extensions  to 
the  basic  serial  system  models.  For  example,  Lambrecht  and  Vander  Eecken 
(1977)  have  extended  the  Zangwill  algorithm  to  include  capacity  constraints 
for  facilities  in  series.  Also, for  the  serial  system  with  concave  costs 
and  storage  costs  non-decreasing  in  order  of  facilities  and  production 
costs  non- increasing  in  time.  Love  (1972)  has  demonstrated  that  an 
optimal  schedule  has  the  property  that  if  in  a given  period,  a facility 
produces,  then  its  successor  does  also.  This  nested  structure  is 
exploited  in  an  algorithm  for  finding  an  optimal  schedule. 
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In  a multi-echelon  assembly  system  each  facility  supplies  at  most 
one  immediate  successor,  but  it  may  be  supplied  by  more  than  one  imme- 
diate predecessor.  External  demands  occur  only  at  the  final  facility. 
Clearly,  the  assembly  of  components  into  subassemblies  and  eventually 
into  one  final  product  is  an  example  of  such  a system.  Also,  the  serial 
system  is  a special  case  of  this  more  general  configuration. 

Crowston,  Wagner  and  Williams  (1973)  investigate  the  continuous 
demand  case  with  instantaneous  production  and  integral  lot  sizes.  Under 
the  assumption  of  a stationary  lot  size  at  each  facility,  it  is  demonstrated 
that  the  optimal  lot  size  at  each  facility  is  an  integer  multiple  of  the 
lot  size  at  the  successor  facility.  This  fact  is  used  in  the  construc- 
tion of  a dynamic  programming  algorithm  for  the  computation  of  optimal 
lot  sizes.  Schwarz  and  Schrage  (1975)  omit  the  assumptions  of  instan- 
taneous production  and  an  integral  lot  size  at  the  final  stage  and 
present  a heuristic  solution  technique  and  a branch  and  bound  algorithm 
to  determine  optimal  lot  sizes.  In  addition,  Schwarz  and  Schrage  assume 
that  no  lot  splitting  is  permitted.  That  is,  no  part  of  a lot  produced 
at  a facility  is  available  for  use  by  its  successor  facility  until  the 
entire  lot  has  been  completed . 

For  the  case  of  dynamic  demand  in  a multi-stage  assembly  system, 
Crowston  and  Wagner  (1973)  present  a dynamic  programming  algorithm 
and  a branch  and  bound  algorithm  to  determine  optimal  lot  sizes  that 
take  advantage  of  Love  type  nested  schedule  for  a particular  cost  structure. 

In  a multi-echelon  arborescent  system,  each  facility  is  supplied  by 
at  most  one  immediate  predecessor,  but  it  may,  in  turn,  supply  more  than 
one  immediate  successor.  External  demands  occur  only  at  the  facilities 
without  successors.  A national  warehouse  which  supplies  a group  of 
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regional  warehouses,  each  of  which,  in  turn,  supplies  a group  of  retailers 
is  an  example  of  such  a system.  An  arborescent  system  with  only  2 
echelons  where  the  upper  echelon  consists  of  only  a single  facility  is 
referred  to  as  a one  warehouse,  N retailer  system. 

A limited  amount  of  research  has  been  conducted  concerning  one 
warehouse,  N retailer  and  general  arborescent  systems.  For  the  dynamic 
demand  case,  Veinott  (1969)  discusses  the  possible  extension  of  Zangwill's 
algorithm  to  problems  with  an  arborescent  structure.  Kalymon  (1972) 
presents  a decomposition  algorithm  which  reduces  the  arborescent  problem 
to  a series  of  single  stage  problems  at  each  lowest  echelon  facility. 

The  dynamic  demand  case  for  the  one  warehouse , N retailer  system  will 
be  investigated  further  in  Chapter  3. 

For  the  continuous  demand  case,  Schwarz  (1973)  proposes  the  use  of 
"singl.i  cycle"  and  "separate  retailing"  policies  for  the  one  warehouse, 

N retailer  system.  Attempts  have  been  made  (Graves  and  Schwarz,  1977) 
to  relate  these  policies  to  more  general  arborescent  structures  and  to 
determine  the  optimality  of  suon  solutions.  A complete  discussion  of 
these  policies  will  be  presented  in  Chapter  2 along  with  a theoretical 
foundation  detailing  their  relation  to  optimality. 


CHAPTER  2 


CONTINUOUS  ONE  WAREHOUSE,  N RETAILER  PROBLEMS 

2.1  Introduction 

In  this  chapter  we  examine  the  continuous  one  warehouse , N retailer 
production  inventory  system.  We  begin  by  introducing  the  basic  assumptions 
as  well  as  some  previously  solved  special  cases.  Schwarz  (1973)  has 
proposed  two  basic  production  policies  called  separate  retailing  and 
single  cycle . We  will  reexamine  both  of  these  policies  with  particular 
emphasis  on  algorithms  to  determine  the  optimal  single  cycle  policy.  We 
will  also  carefully  examine  the  optimality  of  single  cycle  policies  for 
several  wider  classes  of  production  plans. 

2.2  Basic  Model 

The  one  warehouse , N retailer  system  is  a special  case  of  the 
arborescent  production  inventory  system  discussed  in  Chapter  1.  The  lower 
echelon  consists  of  N "retailers",  that  is,  those  facilities  that  directly 
experience  the  external  demand  and  the  upper  echelon  consists  of  one  "ware- 
house", which  is  the  predecessor  facility  to  all  the  retailers.  Raw 
materials  enter  the  warehouse , where  some  processing  takes  place . Units 
flow  from  the  warehouse  to  the  retailers  where  final  processing  takes  place 
and  these  units  then  flow  out  of  the  system  to  satisfy  external  demand. 

We  will  refer  to  the  warehouse  as  facility  0 and  we  will  refer  to  the  N 
retailers  as  facilities  1 through  N respectively  (Figure  2.1).  This 
structure  of  facilities  is  common  to  both  industrial  and  military  situations. 
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EXTERNAL  DEMANDS 

FIGURE  2.1  ONE  WAREHOUSE,  N RETAILER  STRUCTURE 


In  the  Air  Force,  for  example,  the  lower  echelon  facilities  are  referred 
to  as  bases  and  the  upper  echelon  as  a depot  (Muckstadt,  1977). 

In  the  continuous  demand  case,  the  demand  at  retailer  i is  assumed 
to  occur  at  a known  constant  rate  . The  demand  rate  may  differ  from 
one  facility  to  another.  We  will  assume  that  the  warehouse  experiences 
no  external  demand  and  thus  the  demand  at  the  warehouse  occurs  at  a 
constant  rate,  Dg,  where 

D = f D . (2.1) 

0 i=l  1 

At  each  facility  there  is  a set-up  cost  and  a holding  cost.  In 
particular,  let 


= Set-up  cost  incurred  to  produce  a lot  at  facility  i 
and  h^  = Echelon  holding  cost  per  unit  held  at  facility  i, 


*: 
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where  i = 0,1,2,...,N. 

It  is  important  to  point  out  that  the  echelon  holding  cost  is  charged 
against  the  echelon  inventory,  a concept  first  introduced  by  Clark  and  Scarf 
(1960).  For  the  retailers,  the  echelon  inventory  is  the  on-hand  inventory  in 
the  usual  sense.  For  the  warehouse,  the  echelon  inventory  is  equal  to  the 
inventory  held  anywhere  in  the  system.  The  use  of  echelon  inventory  in 
the  model  permits  some  very  convenient  simplifications  but  it  has  been 
pointed  out  that  the  use  of  echelon  inventory  instead  of  the  more  typical 
on-hand  inventory  may  be  confusing  to  anyone  who  wants  to  apply  the 
model  in  practice  ( Szendrovits , 1978).  Schwarz  and  Schrage  (1978)  have 
demonstrated  that  a policy  which  is  optimal  under  the  echelon  stock 
charging  scheme  is  also  optimal  under  a charging  scheme  based  upon  the 
actual  on-hand  inventory. 

With  the  cost  structure  defined  in  the  above  way,  we  can  relax  the 
restriction  that  only  the  retailers  directly  experience  the  external  demand. 
If  the  external  demand  rate  at  the  warehouse  is  D^,  one  can  add  to  the 
problem  an  additional  retailer  with  demand  rate  D^,  a set-up  cost  of  0 
and  an  echelon  holding  cost  of  0.  This  creates  a one -warehouse , N+l 
retailer  problem  with  no  external  demand  at  the  warehouse  which  is  equiva- 
lent to  the  original  problem. 

The  assumption  that  no  lot  splitting  is  permitted  is  common  to  much 
of  the  research  on  multi-stage  production  inventory  systems  (Schwarz 
and  Schrage,  1975).  The  no  lot  splitting  assumption  states  that: 

(1)  No  part  of  a lot  produced  at  a retailer  may  be  used  to  satisfy 
external  demand  until  the  entire  lot  is  produced  and  (2)  No  part  of  a 
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lot  produced  at  the  warehouse  may  be  shipped  to  any  retailer  until  the  entire 
lot  is  produced.  This  assumption  simplifies  the  computational  effort  and 
is  quite  realistic  in  many  cases.  For  example,  the  warehouse  and 
retailers  may  be  physically  far  apart  and  it  is  logical  that  units  will 
be  shipped  in  bulk  rather  than  individually.  For  the  present,  we  will  also 
assume  that  no  lot  splitting  is  permitted,  although  we  will  examine  the 
alternative  in  a later  section. 

In  addition,  we  will  assume  that  the  delivery  times  between  the  ware- 
house and  the  retailers  and  between  the  retailers  and  the  external  customers 
are  constant.  Given  this  assumption  along  with  no  lot  splitting,  we  can 
also  assume  that  production  is  instantaneous.  The  results  presented  can 
easily  be  extended  to  the  case  of  finite  production  rates . Under  no  lot 
splitting,  no  part  of  a lot  can  be  used  until  the  entire  lot  has  completed 
production  and  the  flow  of  units  through  the  system  will  remain  unchanged 
if  we  assume  that  production  is  instantaneous.  Finally,  due  to  the  zero 
delivery  times  and  instantaneous  production  rates,  units  can  be  made 
available  to  satisfy  external  demand  at  time  zero.  Thus,  we  may  also 
assume  that  the  initial  system  inventory  is  zero. 

2.3  Properties  of  an  Optimal  Solution 

We  will  begin  our  investigation  of  production  plans  for  the  one 
warehouse , N retailer  production  inventory  system  by  introducing  some 
properties  that  optimal  policies  must  possess  (Graves  and  Schwarz,  1977): 

PI.  Each  facility  produces  (orders)  only  when  its  on-hand  inventory  is 
zero. 

P2.  The  warehouse  produces  only  when  at  least  one  of  its  successors 
produces . 
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P3.  For  each  retailer,  all  lot  sizes  produced  in  the  time  interval 
between  successive  production  runs  at  the  warehouse  are  equal. 
Properties  P1-P3,  stated  in  a slightly  more  general  fashion,  will 
hold  for  the  optimal  policies  of  general  arborescent  systems  as  well. 

The  proofs  of  properties  P1-P3  may  be  found  in  Schwarz  (1973). 

We  will  define  a simultaneous  production  point  for  a retailer  j 
as  a point  in  time  when  stage  j and  the  warehouse  produce  simultaneously. 
There  is,  clearly,  a simultaneous  production  point  at  time  zero  for  all 
retailers  in  any  feasible  plans  since  we  have  assumed  initial  inventories 
of  0.  In  order  to  ins&re  the  existence  of  a simultaneous  production 
point,  after  time  zero,  it  is  necessary  to  assume  that  h^  > 0.  If 
hQ  = 0,  the  optimal  policy  for  the  warehouse  is  to  produce  an  infinite 
lot  size  at  time  zero,  regardless  of  the  policy  at  any  of  the  retailers. 
Thus , if  hQ  = 0 , the  one  warehouse , N retailer  problem  becomes  simply 
N independent  one-retailer  problems  which  have  been  thoroughly  analyzed 
(Hadley  and  Whitin,  1963).  From  this  point  on,  we  will  assume  that 
hQ  > 0.  Given  this  assumption,  we  can  add  the  following  properties: 

P4.  The  time  interval  between  successive  simultaneous  production  points 
is  finite  for  each  retailer. 

P5 . The  lot  sizes  for  each  retailer  are  non-decreasing  between  its 
successive  simultaneous  production  points. 

Graves  and  Schwarz  (1977)  present  proofs  for  properties  P4  and  P5. 

In  the  proof  for  P4,  they  incorrectly  assume  that  the  greatest  lower 
bound  on  an  infinite  set  of  positive  numbers  is  positive.  A revised 
proof  based  in  part  on  an  incorrect  proof  by  Graves  and  Schwarz  (1978) 
is  presented  below. 
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Proof  of  P4 . Consider  some  retailer  violating  P4  in  an  optimal  policy  P'. 
For  simplicity,  we  will  assume  that  the  only  simultaneous  production  point 
for  retailer  j is  at  time  zero.  (If  this  is  not  the  case,  a simple  shift 
of  the  time  axis  will  make  the  proof  valid.)  Let  I^(t)  be  the  echelon 
inventory  of  stage  j at  time  t.  Let  (t^jt^.t^,. . .)  be  the  times  when 
production  takes  place  at  the  warehouse  in  P'.  (The  existence  of  such  a 
countable  set  of  production  points  is  assured  if  the  set-up  cost  at  the 
warehouse  or  at  all  the  retailers  is  strictly  positive.)  Let 
I1?111  = inf{l.(t.)}.  That  is,  I.(t.)  is  the  amount  of  inventory  at 


I i 


retailer  j at  a point  in  time  when  there  is  production  at  the  warehouse. 


There  are  two  possible  cases:  Case  1,  in  which  I™in  > 0;  and  Case  2,  in 
which  I?ln  = 0.  The  proof  for  Case  1 is  presented  in  Graves  and  Schwarz 
(1977). 

The  proof  for  Case  2 will  construct  a plan  P with  no  greater  cost 
than  the  assumed  optimal  policy  P'  and  with  a simultaneous  production  point 
after  time  0.  In  P,  we  will  reduce  the  amount  produced  at  the  retailer 
in  P'  for  a particular  point  in  time.  This  reduction  will  force  a group 
of  subsequent  productions  to  be  shifted  earlier  in  time.  This  shift  will 
create  a simultaneous  production  point  in  P. 


In 


particular,  define  Q,  = sup{Q.(t.)}  where  t.  are  the  times 


‘j  i' 


when  production  takes  place  at  j and  Q^(t^)  is  the  lot  size  at  retailer 
j at  time  t^  in  plan  P'.  Incidentally,  < ®;  otherwise,  the  cost 
of  P'  would  be  infinite,  assuming  of  course  the  holding  costs  are  posi- 
tive. Let  T = (h^/hQ)  *Qj/Dj . Intuitively,  Q^/D^  may  be  viewed  as  the 
maximum  time  interval  between  successive  production  runs  at  j . 

Define  t,  to  be  the  smallest  t.  > T.  (Such  a t,  exists  if  the 

K 1 ““  K 

holding  cost  at  the  warehouse  is  non- zero.)  Let  e = min  (l.(t.)}  >0. 

i=l,2,...,k  ] 1 
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(If  e=I.(t)=0,  then  there  is  a simultaneous  production  at  t ). 

I * ° 

Consider  policy  P which  is  identical  to  P'  at  the  warehouse 

except  that  at  time  0 the  warehouse  produces  e less  than  in  P' 

and  at  time  t^  the  warehouse  produces  z more.  The  echelon  holding 

costs  at  the  warehouse  on  the  interval  (0,t, ) in  P are  e-t,  -h 

K K U 

less  than  the  echelon  holding  costs  on  the  same  interval  in  P' . (See 
Figure  2.2.) 

We  will  now  show  that  it  is  possible  to  construct  P at  retailer  j 

such  that  the  increase  in  costs  from  P'  is  no  larger  than  the  decrease 

in  costs  at  the  warehouse  evaluated  above.  Let  x^  be  the  earliest  time 

after  t^  when  production  takes  place  at  retailer  j . Also  let 

x'  = min{x! |Q.(t! ) > e}.  That  is,  x'  is  the  first  production  point  for 
m l 1 j i m r r 

which  retailer  j has  lot  size  greater  than  e.  Since  I1?111  = e,  x'  < t„ 
J ° ] m2 

or  Ijttj)  < e contradicting  the  definition  of  e.  We  will  denote  the 

production  points  at  j,  in  the  revised  plan  P,  by  x^.x^,...  . Up 

until  time  x , P is  identical  to  P' . That  is,  • x.  = x!.  and 
m’  ’ii 

P i P 

Q.(x.)  = Q.(x.)  for  1 < i < m,  where  Q.(x.)  is  the  lot  size,  at 
] i 1 i — : i 

retailer  1,  at  time  x.  in  the  revised  plan  P.  Let  x = x'  but 
J l r mm 

P • 

Q.(x  ) = Q.(x  ) - e.  The  conditions  on  m will  insure  that  this  change 
] m j m 

maintains  feasibility  at  x^.  Since  fewer  units  are  produced  at  x^ 
to  satisfy  the  demand  at  retailer  j , the  next  production  point , Tm+1 
must  occur  earlier  than  T^+1  to  maintain  feasibility.  Since  c less 
units  are  produced  at  x^,  Tm+1  - T^+1  “ e/Dj  to  insure  feasibility. 

In  fact,  if  we  set  = ^or  m < i < * then  all  production 

points  x^,  for  m < i <^  *.,  are  shifted  earlier  in  time  than  the 
corresponding  production  point  in  P to  maintain  feasibility  and  they 


16 


are  all  shifted  back  by  the  same  amount,  e/D..  That  is,  t.  = t.  - e/D. 

] ’ll] 

for  m < i <_  i.  Finally,  we  add  e units  to  the  production  at  x^. 

P i < 

That  is,  Q^(x^)  = Qj(x^)  + e*  On  the  interval  P i-s  identical 

to  P' . 

To  clarify  this  construction,  at  retailer  j,  consider  Figure  2.2. 

In  this  particular  case,  it  is  assumed  that  m=4,  k=4,  1=7  and 

e = Ij(t2>.  The  holding  costs  under  both  P1  and  P are  illustrated, 

with  the  holding  costs  under  the  revised  plan  P denoted  by  the  shaded 

area.  P is  identical  to  P'  on  (0,x^)  and  on  (x^,®).  At  x^,  e 

fewer  units  are  produced  in  P and  thus,  retailer  j reaches  zero 

inventory  at  x = x'  - e/D..  x and  x are  also  shifted  e/D.  units 
o b 3 b / 3 

» | ... 

back  in  time  from  xc  and  x_.  At  x_,  e additional  units  are  produced 

b II 

and  thus  at  x^  plan  P is  again  in  phase  with  P1. 

Let  us  denote  t ..  as  the  production  point  in  Pf  at  the  warehouse 

for  which  = I.(t  . ) = e.  At  (t  . + e/D.),  in  P' , retailer  j 

] ] min  min  ] 

has  0 inventory  and  thus,  by  the  feasibility  of  P' , (t  . + e/D.)  = x! 

min  ] l 

for  some  m < i < i,.  Therefore,  t . = x!  - t.  J.  = x.  by  the  construe - 

— min  i ] l 

tion  of  P.  Thus,  t . = x.  is  a simultaneous  produce >>  point  after 

min  l 

time  0 in  the  plan  P.  In  Figure  2.2,  I^(t2)  = e and  t^  = x^  is 

a simultaneous  production  point  in  P. 

It  is  obvious  from  Figure  2.2  that  the  holding  cost  in  P is  less 
than  the  holding  cost  in  P1  on  the  interval  (0,x^)  and  equal  to  the 
costs  in  P*  on  (x^,®).  In  general,  the  only  increase  in  costs,  from 
P*  to  P,  is  on  the  interval  p»  QjCx^)  additional 

units  are  held  over  the  interval  (x,x^).  Therefore  the  increase  in  costs 


can  be  written  as 
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h. *(t  -t  )*Q. Ct')  = h.-(e/D  )-Q  (t'). 

JXXJX,  j Jj^ 

Furthermore,  by  definition  Q.(x')  <_Q.,  and  t > T.  Substituting  in 

1 * 3 * 

the  definition  of  T and  rearranging  terms , we  see 

e*t  *h  1 e-h.-Q./D  > h. • (e/D. ) *Q ,(t' ) . 

KU  J J J J j I ^ 

Therefore,  the  decreases  in  holding  cost  from  P'  to  P exceed  the 
increases  in  holding  cost  and  since  P has  the  same  number  of  set-ups  as 
P' , P is  no  more  costly  than  P'.  By  repeated  use  of  this  construction 
plan,  we  can  show  that  there  is  a plan,  at  worst  identical  in  cost  to 
the  optimal  plan,  for  which  the  time  interval  between  successive  simul- 
taneous production  points  is  finite.  □ 

2.4  Renewal,  Single  Cycle  and  Stationary  Policies 

We  define  a renewal  policy  to  be  one  in  which  there  exist  points  in 
time  t,  0 < t < 00 , at  which  all  stages  produce  simultaneously.  If 
there  exists  a renewal  policy  with  a point  at  which  all  stages  produce  and 
some  stage  j has  non-zero  inventory,  e,  we  can  reduce  the  costs  of 
this  plan  by  reducing  the  previous  production  lot  at  stage  j by  e . 

Thus,  we  will  always  assume  that  at  each  renewal  point,  all  stages  have 
zero  inventory.  Given  this  assumption,  the  system  at  a renewal  point 
appears  exactly  as  if  we  were  starting  the  system  again  at  time  zero. 

If  a renewal  policy  is  the  overall  optimal  policy  and  the  first  renewal 
occurs  at  T,  we  will  clearly  have  renewals  every  T time  units  and 
the  policy  will  be  the  same  between  each  pair  of  renewals.  Furthermore, 
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in  the  best  renewal  policy  the  retailers  produce  the  same  lot  size  each 
time  they  produce.  Alternatively,  we  can  state, 

Lemma  2.1.  The  optimal  renewal  policy  is  stationary  at  the  retailers. 
Proof  of  Lemma  2.1. 

Suppose  the  optimal  renewal  policy  is  not  stationary  at  retailer  j 
i.e.  between  a pair  of  renewals,  retailer  j produces  consecutive  lots 
of  sizes  and  where  i-  Q^.  The  holding  cost  in  this  policy 

is : 

k(vv  * kcvv  • kkkv 

Look  at  the  revised  policy  where  we  produce  two  lots  of  size  (Q^+Q^)/^. 
This  policy  certainly  is  feasible  if  the  original  policy  is  feasible. 

The  cost  of  this  policy  is: 

2(j<((J1tQ2)/2)(<(Q1t<32)/2)/D1))  = ((Q1+Q2)/2)2/D1. 

/VQA  2 Q2«)2 

So  if  I — - — J < — - — then  the  revised  policy  is  cheaper  than  the 
original.  This  is  clearly  the  case.  □ 

In  a renewal  policy,  the  warehouse  produces  at  every  renewal  point 
by  definition.  Furthermore,  suppose  that  at  every  warehouse  production 
point,  we  also  have  a renewal.  This  special  type  of  renewal  policy  is 
referred  to  as  a single  cycle  policy.  The  best  single  cycle  policy  is 
clearly  stationary  at  the  retailers  and  at  the  warehouse. 


20 


(2.2) 

(2.3) 

(2.4) 

(2.5) 

where  n.,  j > 1,  is  defined  as  the  integer  number  of  lots  delivered  to 
, stage  j from  the  warehouse  and  T is  the  length  of  time  between  renewals. 

The  equation  (2.2)  is  the  expression  for  the  average  cost  per  unit 
time  for  the  single  cycle  policy  (n^^ , . . . ,nN) . Clearly,  if  the 
warehouse  sets  up  every  T time  units,  the  average  set  up  cost  per  unit 
time  is  K /T.  Retailer  i sets-up  n.  times  every  T time  units  so 

I 1 

the  average  set-up  cost  per  unit  time  for  retailer  i is  n.K./T.  If 

I 

retailer  i produces  a lot  of  size  Q.  at  each  production,  then  because 

1 Q. 

1 

of  the  constant  demand  rate  D.,  the  retailer  holds  units  on  the 

f 1 ^ 

average.  The  total  demand  at  retailer  i during  the  cycle  length  T 

T-D. 

is  T-D..  Since  n.  equal  lots  are  produced  by  Lemma  2.1,  Q.  = — . 

; i 

* The  average  holding  cost  per  unit  time  at  retailer  i can  therefore 

I 

be  written  as  h^T-D^/2n^.  The  warehouse  holding  cost  can  be  expressed 
siu»i1arly  because  the  echelon  inventory  is  drawn  down  at  the  constant 
rate,  DQ.  The  validity  of  using  the  average  cost  per  unit  time  formu- 

k 

lation  has  been  verified  by  Schwarz  (1973). 


In  Figure  2.3,  a two  retailer  single  cycle  policy  is  illustrated 

where  n.  =2  and  n„  = 3. 
i 


«!■» 


The  optimal  single  cycle  policy  may  be  found  by  solving: 

N /n.K. 

Minimize  z = \ + 

j=0\  T 

subject  to  T >_  0 , 


n.  >_  1,  integer  for  j = 1,...,N, 
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Let  us  consider  the  implications  of  relaxing  our  assumption  that 
no  lot  splitting  is  permitted.  In  this  case,  we  can  also  assume  a finite 
production  rate  p^  at  each  facility  i = 0,1,..., N.  If  we  start  with 
zero  initial  inventories,  to  insure  feasibility,  we  must  have 


pQ  > Pi  t ...  t PN  and 

p.  > D.  for  all  facilities  i. 
1—1 


In  this  case,  rather  than  having  an  immediate  jump  to  a height  of 

Ql , as  in  the  no  lot  splitting  case , the  echelon  inventory  at  facility 

i increases  for  Q./p.  time  units  at  a rate  of  (p.-D.)  and  then 
1 1 11 

decreases  at  rate  D.  for  Q./D.  - Q./p.  time  units.  If  the  cycle 
length  is  T,  then  Qi  = TD.j/n£  for  the  single  cycle  policy 
(°1  ,n2  * ' ‘ ’ ,njj>  * 't^^s  oase , the  average  echelon  inventory  held  per 

unit  time  at  facility  i can  be  expressed  as 

TD.  D. 


(1-P.) 


where 


1 


Thus,  the  optimal  single  cycle  policy,  with  lot  splitting  allowed,  may 
be  found  by  solving 


r 


i 
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N /n.K.  h.TD.  > 

Minimize  z ~ l (^~T  + "2'n.^  (1_Pjy 


subject  to  T > 0, 


n0  = li 


n^  >_  1,  integer  for  j = 


and  p.  = — 
1 Pi 


Since  p^,  is  a constant  we  can  define 


hj  ■ 


and  the  optimization  problem  above  reduces  to  precisely  the  one  described 
by  (2-2)  — (2.5),  that  is,  finding  the  optimal  single  cycle  policy  with  no 
lot  splitting  permitted. 

The  following  results  concerning  single  cycle  policies  have  been 

established  (Schwarz,  1973). 

•? 

Theorem  2.1.  For  the  one  warehouse,  one  retailer  system,  the  optimal 
policy  is  a single  cycle  policy.  □ 


Theorem  2.2.  For  the  one  warehouse,  N retailer  system,  where  the  N 
retailers  have  identical  costa  and  demand  rates,  the  optimal  policy  is 
a single  cycle  policy.  0 


We  can  derive  more  general  relationships  concerning  the  optimal  single 


cycle  policies  and  the  retailer  costs. 
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Theorem  2.3.  If  h.D./K.  < h.D./K.,  then  n*  < n*  in  the  optimal  single 
i i i 3 3 3 i-l  * 

cycle  policy. 

Proof  of  Theorem  2.3. 

Suppose  we  have  the  optimal  single  cycle  solution  (n*,n*, . . . ,n*) 

and  an  optimal  cycle  length  T*  such  that  l^D  /K.^  < h2D2^K2  but 

n“  > n*.  The  cost  of  this  policy  is  less  than  or  equal  to  that  for  all 

policies  with  the  same  cycle  length  but  a different  vector  of  n^'s.  In 

particular,  the  cost  of  the  optimal  policy  is  smaller  than  or  equal  to 

the  cost  of  the  solutions  (n*,n*,n*  . . . ,n*)  and  (n* ,n* ,n*, . . . ,n*) 

1 I o N 2 2 o N 

both  with  cycle  length  T*.  If  we  were  to  write  out  the  cost  function 
(2.2)  for  these  3 plans  with  cycle  length  T*,  the  only  terms  that 
differ  involve  retailers  1 and  2.  Thus 


and 


niKl  n2K2  T*hiDi  T*h2D2 

n*K1  + n*K2 

T*(h^D^  + h2D2) 

14  ' T*  ' 2n*  ' 2n* 

l 

— 

* 

"X  “X  T*hlDl  T*h2D2 

n*K,  + n*K_ 

_ 2 1 2 2 

T*(h1D1  + h2D2) 

T*  ' T*  ' 2n*  ' 2n* 

— t* 

* 

(2.6) 


(2.7) 


Thus 


and 


n*K2  T*h2D2  n*K2  T*h2D2 


n*K  T*h  D n*K  T*h  D 
1*1  11  < 21  11 

T*  „_ft  - T*  + „_ft 


2n, 


2n„ 


Thus 


ft 

*\ 

n,  - 

■n„  \ 

1 

_2  ] 

ft 

ft  J 

n2ni  / 

1 (n.-n.) 


2 
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and 


K T*h  D / n*-n"^ 

(n*-n*)-i  < L-Lf  1 2 

^ 1 2;T*  - 2 \ * * 


Dividing  both  sides  of  both  inequalities  by  n”  - n*  > 0 we  get: 


h2°2  2n2ni  hlDl 
K2  “ (T*)2  “ K1 


h2D2  hlDl 

This  contradicts  our  original  assumption  that  — — > — — . □ 

K2  K1 

Theorem  2.4.  If  h.D./K.  = h.D./K.,  then  n*  = n*  in  an  ODtimal 
l l i 3 ] : i 1 H 

single  cycle  policy. 


Proof  of  Theorem  2.4. 

Suppose  we  have  the  optimal  single  cycle  solution  (n*,n*, . . . ,n*) 

12  N 

and  an  optimal  cycle  length  T*  such  that  h1D1/K1  = h2°2/K2  but 
nl  ^ n2‘  cost  'this  policy  is  less  than  or  equal  to  that  for  all 

policies  with  the  same  cycle  length  but  a different  vector  of  n^'s  and 
as  in  the  proof  of  Theorem  2.3  we  again  arrive  at  equations  (2.6)  and 
(2.7).  Thus 


and 


n2%  T*h2D2 

< 

T‘ 

_ t* 

T*hlDl 

< 

~fr'  ft 

2n* 

~ T* 

T*h.D 
2 2 

2n * 


T*h^D^ 

2n* 


(2.8) 


(2.9) 


Rewriting  (2.8)  we  have 
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Since  by  assumption 


By  (2.9) 


(n*-nj)  < (T *Y 


h2D2/K2  = hlDl/K 


(n*-n*)  1 (T*)2' 


(2.10)  becomes 


(n*-n*) , 


(2.10) 


(2.11) 


(2.12) 


So  clearly  for  (2.11)  and  (2.12)  to  hold  simultaneously,  both  relations 

must  be  equalities  which  implies  that  (2.8)  and  (2.9)  are  also  equalities. 

> 

Then  clearly  we  can  change  the  optimal  policy  from  (n*,n*,n*,. . . ,n*) 

12  3 N 

to  either  (n*,n*,n*,. . . ,n*)  or  (n*,n*,n*, . . . ,n*)  without  any  increase 
in  cost.  □ 


2.5  Optimal  Cycle  Length  and  Separate  Retailing 

We  will  begin  by  examining  the  cost  for  a single  cycle  policy  in 
greater  detail.  Recall  (2.2)-(2.5) 


N f n.K.  h.TD.X 


T > 0, 


n0  = 1* 


nj  i1*  j = I,-  • • »N. 
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For  fixed  values  of  ru , j = 1,...,N,  the  optimal  T,  denoted  T*(n^), 
satisfies : 


T*(n . ) = 
J 


2 J"  n.K. 

j = 0 3 ^ 

N 

l h,D./n. 

L j=0  ^ 3 


1/2 


(2.13) 


and  the  optimal  average  cost  per  unit  time  denoted  C*(n^),  satisfies: 


N N 1/2 

C*(n . ) = [ 2 ( l n.K. )(  [ h.D./n.)] 


: 


j =o  ^ y j =o  J j y 


(2.14) 


Since  (2.2)  is  separable  and  convex  in  the  n^ , we  may  determine  the 
optimal  n_. , j = 1,...,N,  for  fixed  T,  denoted  n\‘(T),  as  the  smallest 
n satisfying 


n.(n.+l)  > T (h.D./2K.). 
d 3 - : : d 


(2.15) 


Given  the  value  C*(n^ ) as  in  (2.14),  it  is  possible  to  show  that  if  two 
problems  have  the  same  basic  cost  structure,  then  the  optimal  solution 
vector  of  n^'s  will  be  the  same.  More  specifically,  we  can  state 

Theorem  2.5.  If  we  have  2 different  one  warehouse,  N retailer  problems 
with  costs  K.,h.,D.  and  K!,h!,D'.,  i = 0,...,N  where 

l l i i l l’ 


K.  = oK.  i = 0 , . . . ,N 
i i 


and 


i 


hi'Di  = 0,VDi  L = 
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J.  m 

then  the  optimal  single  cycle  solution  vectors  (n^,...,n*)  and 

(n*' ,. . . ,n*' ) are  identical. 

1 N 

Proof  of  Theorem  2.5. 

Suppose  (n*,...,n”)  is  the  optimal  solution  for  the  first  problem. 

Then 

C(n*,...,n*)  < C(n_ , . . . ,n„)  for  n.  = 1 and  for  all 
1 N — 1 N 0 

n1,...,nN,  integer,  ■>  1 (2.16) 

and  by  (2.14) : 

[2(  X n*K.)(  f h.D./n*)]1/2  < [2(  f n.K.)(f  h.D./n.)]1/2 

j = 0 3 ^ j=0  J j=Q  ] i j=Q  J 1 3 

for  nQ  = 1 and  for  all 

n. ,...,n„,  integer,  > 1.  (2.17) 

IN  — 

1/2 

Multiply  both  sides  of  (2.17)  by  (a8)  . Then  regrouping  terms  yields 

N N . N N ... 

C 2 ( l n*aK.)(  l Bh.D./n?)]  ' < [2(  £ n.aK.X  l Bh.D./n. )]  ' 

j = 0 3 J j=0  3 J j =0  3 : j-0  3 ] 3 

for  nQ  = 1 and  for  all 
n1,...,nN,  integer,  >.  1. 

This  clearly  demonstrates  that  (n*,...,n*)  is  also  the  optimal  solution 
for  the  problem  with  costs  K^,h|,D^. 

These  two  problems  do  not  have  identical  ratios  but  the  relationships 
between  these  ratios  within  each  problem  are  the  same.  0 
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We  will  now  examine  the  one  warehouse,  one  retailer  problem  more  com-  < 

pletely.  For  one  retailer,  i.e.  N = 1,  DQ  = D and  (2.14)  becomes 

C*^)  = [2(KQ+n1K1)(h0D1+h1D1/n1)]1/2.  (2.18) 

Given  the  convexity  of  the  expression  within  the  brackets  of  (2.11),  the 

JL 

optimal  n^,  denoted  n^,  is  chosen  as  the  smallest  integer  satisfying: 

n(n+l)  >_  KQh  /K  h . (2.19) 

Thus  we  can  explicitly  solve  the  one  warehouse,  one  retailer  case 
to  find  the  best  single  cycle  policy,  which  by  Theorem  2.1  is  the  overall 
optimal  policy. 

This  suggests  an  alternative  policy  for  the  one  warehouse,  N retailer 
system.  If  we  split  the  one  warehouse,  N retailer  system  into  N one 
warehouse,  one  retailer  systems  and  solve  each  of  these  one  retailer  systems 
separately,  we  will  arrive  at  the  policy  called  separate  retailing. 

It  is  clear  that  the  separate  retailing  policy  is  the  overall  optimal 
solution  for  the  special  case  where  = 0.  If  = 0,  then  the 
on-hand  inventory  at  the  warehouse  is  always  zero  in  the  overall  optimal 
solution.  Therefore,  the  warehouse  echelon  holding  cost,  hQ,  is  only 
charged  against  the  inventory  held  at  the  retailers.  In  this  case,  we 
can  solve  N independent  single  facility  problems  with  the  i**1  such 
problem  having  a set-up  cost  of  and  a holding  cost  of  (h^+hQ), 

and  then  schedule  set-ups  at  the  warehouse  whenever  they  are  necessary. 

By  setting  KQ  to  zero  in  (2.18)  and  (2.19)  it  is  clear  that  the 
procedure  outlined  above  is  simply  separate  retailing. 
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2.6  An  Algorithm  to  Find  the  Optimal  Single  Cycle  Polic 


In  this  section,  we  will  present  a method  that  is  guaranteed  to  find 
the  optimal  single  cycle  policy.  Goyal  (1974)  has  presented  an  algorithm 
that  determines  the  optimum  packaging  frequency  of  items  jointly 
replenished.  A number  of  items  are  manufactured  jointly  and  then 
packaged  with  different  frequencies.  Goyal' s method  determines  an 
N,  the  frequency  of  manufacturing  set-ups  for  the  product,  and  a , 
the  frequency  of  manufacturing  set-ups  for  the  product/ frequency  of 
packaging  set-ups,  to  minimize 


11  . ii 

C(N,K.)  = N(S  + l S./K  ) + l h.-Q.K 
3 j =0  3 3 2N  j = 0 3 3 3 


(2.20) 


where 


n+1  = number  of  items , 


S = manufacturing  set-up  cost, 

Q_.  = demand  per  year, 

hj  = stock  holding  cost  per  unit  per  year , 
Sj  = cost  of  a packaging  set-up. 


It  should  be  noted  that  N is  continuous  and  is  discrete.  If  we 

redefine , 


N = T,  S = 0,  Sj  = h^D V2,  K.  = ru  , hjQj/2  = Kj  and  n = N 


then  (2.20)  becomes 
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N N 

C(T,n. ) = T(  l h D /2n  ) + k £ K.-n.) 

J j 3 J 3 1 j_Q  J J 

N 

= J (n.K./T  + h.D.T/2n.).  (2.21) 

j=0  33  33  3 

This  is  simply  the  cost  of  the  single  cycle  policy  given  in  (2.2)  provided 
nQ  = 1.  Thus,  the  algorithm  presented  by  Goyal,  with  suitable  redefini- 
tion of  terms,  will  find  the  optimal  single  cycle  policy. 

In  order  to  understand  Goyal' s method,  we  will  recall  from  (2.15) 

A 

that  given  T,  the  best  n_. , denoted  n^.(T),  has  the  property  that 

T2h.D. 

n*(T)(n?(T)  - 1)  < — < n*(T)(n*(T)  + 1)  (2.22) 

1 1 • — * 1 1 

1 

Thus , 

ni(n*(T)(n*(T)  - 1))1/2  < T <_  ni(n*(T)(nJ(T)  + 1))1/2  (2.23) 


where 


(2K./h.D. )1/2. 
1 11 


(2.24) 


For  each  retailer  i,  we  can  summarize  this  information  as  in  Table  2.1 

which  gives  the  range  of  the  cycle  length  T for  which  n^(T)  remains 

constant.  We  can  put  all  these  tables  together  to  form  one  master  list 

of  intervals,  I ,1  ,1  ,...  where  the  values  of  n.  remain  unchanged  in 
1 2.  *3  3 

each  interval.  We  can  then  compute  the  cost  of  each  permissible  solution 
by  going  through  this  master  list  interval  by  interval  each  time  com- 
paring the  cost  to  the  best  value  found  so  far.  Since  at  each  step  we 
only  change  a single  n.  by  one,  we  will  eventually  hit  upon  the  optimal 


solution. 


1 
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T 

n . (T)  = 1 

l 

n.  (T)  = 2 

l 

n . (T)  = 3 
l 

n.  (T)  = 4 
l 

Lower  bound 

0 

ni(2)l/2 

n.(6)1/2 

^(20, 1/2 

Upper  bound 

„.<2>1/2 

n.(6)1/2 

„.C12)1''2 

ni(30)1/2 

Lower  and  Upper  Bounds  for  T for  Different  Values  of  n. ($") 


Table  2.1 


It  is  clear  from  (2.13)  that  T*(n.)  increases  as  the  n.  increase. 

3 1 


In  particular. 


T(l,l,...,l)  < TCn^n^,  • . . , ) 


if  at  least  one  n.  > 1 and  all  the  other  n.  >1.  Thus  we  can  begin 
i 3 - 


the  Goyal  procedure  with  cycle  length  T(1 ,1, . . . ,1) . To  halt  the 
procedure  we  will  use  an  upper  bound  on  the  value  of  T*.  Clearly, 


C(n*,. . 


N n.K.  T*h.D. 

‘ li 


K T*hADA 

_ _0  0 0 . t i l , 

• >nN)  T*  + 2 if1  ~W  2n* 


KA  T*hADA  N . /A 

***-£**&  (2KihiDi> 


(2.25) 


This  is  the  case  because  if  we  relax  the  integer  restriction  on  the  n^, 


then 


,1/2 


n.(T)  = T(hiDi/2Kj,)  ' = T/ni 


(2.26) 
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where  is  defined  as  in  (2.24),  and  we  can  substitute  this  value  for 

1/2 

n^(T)  into  the  cost  expression  and  get  (2Kjh^D^)  . In  addition. 


so 


C(l,l,...,l)  > C(n*  ...,n*) 

— 1 N 


C(l,l,...,l)  > Kq/T*  + T*hQD0/2  + l (2K.h.D.) 

i=l 


1/2 


(2.27) 


Define 


N 


A(l,l,l,...,l)  = C(l,l,...,l)  - l (2K.h.D. )1/2. 

i=l  1 1 1 


Thus 


Kq/T*  + T*hQD0/2  < A( 1,1, . . . ,1) . 


The  left  hand  side  of  (2.28)  is  convex  in  T*  and  therefore 


(2.28) 


(2.29) 


T*  < T( 1,1, . . . ,1) 


where 

A(l,l,...,l)  + [ ( A(1 ,1, . . . ,1) )2  - 2K  h D 11/2 
T(l,l, . . . ,1)  = r-g : .(2.30) 

nouo 


This  bound  on  T*  can  be  used  to  halt  the  search  procedure  on  the  n^'s. 

It  should  be  pointed  out  that  Silver  (1975)  has  presented  a simple 
graphical  procedure  to  solve  the  problem  presented  by  Goyal.  This  same 
procedure  is  not  applicable  to  the  single  cycle  problem  because  of  the 
restriction  that  the  warehouse  only  produces  once  in  each  cycle 
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2.7  A Heuristic  Algorithm  to  Find  Good  Single  Cvcle  Policies 


The  Goyal  algorithm,  presented  in  the  preceding  section,  increases 
only  one  n^  by  1 at  each  iteration  and  at  each  iteration  a comparison 
is  made  to  the  incumbent  solution.  In  the  heuristic  method,  we  will 
reduce  the  number  of  comparisons  to  the  incumbent  by  allowing  the  n^ 
to  increase  more  quickly  if  possible.  Given  a set  of  n^'s  we  can 
compute  the  best  cycle  length  for  those  n^'s,  denoted  T*(n^),  using 
the  formula  in  (2.13)  and  given  the  cycle  length  T,  we  can  compute 
the  best  unrestricted  n^(T)  using  (2.26).  The  heuristic  is  essentially 
an  iterative  procedure  that  alternates  between  computing  the  T*(n^)  and 
the  n^(T).  Since  the  minimum  cost  given  a vector  of  n^'s,  denoted 
C*(n^)  and  computed  as  in  (2.14), is  not  a convex  function,  we  could 
converge  upon  a local  minimum  by  simply  iterating  back  and  forth  between 
T*(n^)  and  n^(T).  We  will  avoid  this  problem. by  insuring  that  the 
cycle  length  is  increased  at  each  iteration  which  may  force  us  to  move 
in  a direction  that  does  not  decrease  the  cost,  C*(n^). 

As  in  the  Goyal  method,  we  will  begin  by  evaluating  T(l,l,...,l) 
and  C(l,l,...,l).  We  can  compute  the  upper  bound,  T(l,l,...,l),  on 
T*  using  (2.30)  as  in  the  Goyal  algorithm.  We  will  use  T(l,l,...,l) 

rhiDii 1/2 

to  compute  the  best  unrestricted  n. (T(l,l, . . . ,1) ) = T(l,l,...,l) 

1 2K . 


for  all  i.  We  will  examine  all  the  lattice  points  surrounding  this 
possibly  non- integral  vector  of  n^’s.  By  examining  all  lattice  points 
surrounding  the  unrestricted  n^'s,  the  - n^(T)  as  defined  by  (2.15) 
are  included.  We  will,  however,  only  compute  the  cost  for  a vector  of 
integer  n/s  provided 
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f 


i 


f 


n.  >1  i = 1, . . . ,N, 

l — 


and 


n . > 1 for  some  1 < n . < N , 

3 - ] ~ 


n.  < n.  if  h.D./K.  < h.D./K., 
i - : ill  ] ] ]’ 


n.  = n.  if  h.D./K.  = h.D./K.. 

i 3 ill  333 


Clearly,  we  need  only  check  solutions  that  satisfy  Theorems  2.3  and  2.4. 
All  the  n^  _>  1 from  the  original  statement  of  the  problem  and  there 
is  no  point  in  computing  the  cost  of  (1,1,. ..,1)  again. 

Let 


n^  = smallest  integer  greater  than  or  equal  to  n^(T( 1 ,1, . . . ,1) ) . 

We  can  compute  Kn^n^,...^)  and  check  if  this  cycle  length  exceeds 
the  upper  bound  on  T*  determined  using  (2.30).  If  it  does  not  exceed 
the  upper  bound,  we  will  compute  the  best  unrestricted  n^'s  based  on 
this  larger  cycle  length  and  repeat  the  procedure  outlined  above.  If 
it  does  exceed  the  upper  bound , we  will  halt  the  procedure . It  is  also 
possible  that  when  we  generate  the  new  unrestricted  n^'s  that  all  the 
surrounding  lattice  points  have  already  been  examined.  In  this  case, 
we  will  halt  the  procedure  as  well.  To  improve  the  bounding  procedure, 
we  can  recompute  the  upper  bound  on  T*  whenever  we  find  a better 
incumbent  solution  by  computing  T* ( Incumbent ) using  (2.30). 

To  demonstrate  the  heuristic  procedure , let  us  examine  the  following 
example : 
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Kx  = 99.9,  h1  = 99,  D1  = 1, 

K2  = 99.9,  h2  = 199,  and  D2  = 1. 


(2.31) 


To  begin,  we  will  compute  1/n^  for  all  i,  where  n.  defined  as  in 

1/2 

(2.24).  That  is,  1/n^  = [h^D^/2K^]  . From  Theorem  2.3,  we  know  if 

ft  JU 

1/ n ^ < l/n2J  then  n.^  ± n2‘  ? or  this  example  l/n1  = .7039  and 
l/n2  = *9980  so  we  need  only  examine  those  lattice  points  for  which 
n^  <_  n2*  Next  we  compute 


C(l,l)  = 346.323,  T(l,l)  = 1.15  and  T(l,l)  = 6.266. 


Next  we  compute  the  best  unrestricted  (n^,n2)  for  T(l,l): 

n1(T(l,l))  = l/^-Td.l)  = ,8126  and  n2(T(l,l))  = l/n2-T(l,l)  = 1.15. 

The  only  lattice  point  to  be  examined  is  (1,2)  since  (1,1)  has  already 
been  checked  and 


C( 1,2 ) = 346.727,  T(l,2)  = 1.729. 

Since  T(l,2)  < T(l,l),  we  continue  to  the  next  iteration.  We  compute  the 
best  unrestricted  (n^nj)  for  T(l,2): 


n1(T(l,2))  = 1.21  and  n2(T(l,2))  = 1.72. 
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The  lattice  points  surrounding  (1.21,  1.72)  are  (1,1),  (1,2),  (2,1)  and 

(2.2) .  The  first  two  of  these  four  lattice  points  have  already  been 
examined  and  the  third  can  be  eliminated  since  n > n^  and  by 
Theorem  2.3  this  is  impossible. 

The  only  lattice  point  to  be  examined  is  (2,2)  and 

C( 2 ,2 ) = 347.432,  T(2,2)  = 2.30  < T(l,l). 

We  compute  the  best  unrestricted  (n^.n^)  for  T(2,2): 

n^(T(2 ,2 ) ) = 1.61  and  n2(T(2,2))  = 2.29. 

The  lattice  points  surrounding  (1.61,  2.29)  are  (1,2),  (2,2),  (1,3)  and 

(2.3) .  Only  the  last  two  lattice  points  have  not  yet  been  examined. 
Therefore  we  compute 

C(l,3)  = 365  and  C(2,3)  = 343.13125. 

Clearly  C(2,3)  < C(l,l)  and  (2,3)  becomes  the  new  incumbent  solution. 
We  can  compute  T(2,3)  < T(l,l)  using  (2.30)  and 

T( 2 ,3 ) = 3.056. 

To  insure  that  (n^,n2)  increase  the  most  at  the  next  iteration  we  next 
compute  T( 2 ,3 ) since  (2,3)  is  the  largest  lattice  point  neighboring  the 
non- integer  n/s  from  the  previous  iteration.  We  therefore  compute 

T(2,3)  = 2.912  < T(2,3) 

l 

J 
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and  we  continue  the  algorithm  by  computing 

n (T(2,3))  = 2.05  and  n (T(2,3))  = 2.91. 

The  only  previously  unexamined  lattice  points  are  (3,2)  and  (3,3).  (3,2) 

can  be  eliminated  by  Theorem  2.3  and  so  we  need  only  compute 

C( 3 ,3 ) = 346.84  and  T(3,3)  = 3.457. 

T(3,3)  > T(2,3)  and  so  we  can  halt  the  procedure.  The  best  solution 
found  by  the  heuristic  is  (2,3)  and,  in  fact,  this  is  the  optimal  single 
cycle  policy  for  this  problem. 

It  should  be  pointed  out  that  although  the  example  presented  above 
has  a local  minimum  at  (1,1),  the  heuristic  procedure  does  not  get 
"caught"  by  it.  This  is  due  to  the  fact  that  the  heuristic  does  not 
have  to  move  in  a direction  that  decreases  cost.  Regardless  of  the 
costs,  the  heuristic  increases  the  cycle  length  at  each  iteration. 

2. 8 Computational  Comparisons 

The  Goyal  and  heuristic  algorithms  have  both  been  programmed  in 
order  to  assess  the  quality  of  the  heuristic  and  possible  differences 
in  computational  effort.  The  problems  ranged  in  size  from  2 to  7 
retailers , the  holding  and  set-up  costs  were  uniformly  distributed  on 
the  integers  from  1 to  100  and  the  demand  rates  were  uniformly  distributed 
on  the  integers  from  1 to  10.  The  measure  of  effectiveness  is  the  number 
of  times  the  algorithm  computes  the  cost  of  a complete  plan  and  compares 
that  cost  to  the  value  of  the  incumbent  solution.  To  keep  the  two 
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procedures  on  an  equal  footing,  the  updating  of  the  upper  bound  on  T* 
discussed  for  the  heuristic  is  not  used. 

On  the  average , the  heuristic  procedure  performed  better  than  the 
Goyal  algorithm.  For  the  140  test  problems,  the  Goyal  algorithm 
averaged  8.32  comparisons  to  the  incumbent  the  heuristic  averaged  6.61 
comparisons  per  problem.  More  importantly,  the  heuristic  found  the 
optimal  solution  in  all  140  problems.  There  is,  however,  one  disturbing 
factor.  In  22%  of  the  problems,  the  heuristic  procedure  required  more 
comparisons  than  the  Goyal  algorithm.  In  the  worst  case,  the  Goyal 
algorithm  required  119  comparisons  and  the  heuristic  150.  This  might 
seem  contradictory,  since  the  heuristic  was  constructed  in  such  a way 
as  to  allow  the  vector  of  n^'s  to  increase  more  quickly  at  each 
iteration  than  in  the  Goyal  algorithm.  The  problem  with  the  heuristic 
lies  in  the  examination  of  all  the  lattice  points  surrounding  the  vector 
of  non-integer  n^'s.  For  certain  problems,  this  can  lead  to  a large 
number  of  unnecessary  comparisons  when  using  the  heuristic  procedure. 

As  the  number  of  retailers  increases,  the  number  of  lattice  points  to  be 
examined  at  each  iteration  in  the  heuristic  increases  exponentially. 

The  performance  of  the  heuristic  therefore  declines  in  relation  to  the 
Goyal  algorithm  for  larger  problems.  Fortunately,  a simple  selection 
rule  can  improve  the  computational  efficiency  of  the  heuristic  signifi- 
cantly without  greatly  hurting  its  performance  in  terms  of  finding  the 
optimal  solution.  This  improved  heuristic  will  be  discussed  in  the  next 


section. 
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2.9  An  Improved  Heuristic  Algorithm  to  Find  Good  Single  Cycle  Policies 

In  this  section,  we  will  discuss  a method  to  revise  the  heuristic 
algorithm  presented  in  Section  2.7.  Essentially,  instead  of  examining 
all  the  lattice  points  surrounding  the  non-integer  n^  vector,  we  will 
examine  only  two  particular  lattice  points.  Thus  the  maximum  number  of 
solutions  to  be  evaluated  and  compared  to  the  incumbent  at  each  iteration 
is  two,  regardless  of  the  number  of  retailers  in  the  problem. 

We  would  like  to  ensure  that  the  revised  version  of  the  heuristic 
examines  fewer  solutions  than  the  original  heuristic.  Since  the  maximum 
number  of  solutions  examined  at  each  iteration  in  the  revised  heuristic 

t 

is  fewer  than  in  the  original  approach,  we  need  only  guarantee  that  the 
number  of  iterations  in  the  revised  heuristic  does  not  increase  from 
the  original.  Recall  that  in  the  original  heuristic,  the  largest  lattice 
point  neighboring  the  non-integer  vector  of  n^'s  is  'the  one  which  is 
used  to  get  the  next  value  of  T,  the  cycle  length.  This  same  rule  will 
be  used  in  the  revised  heuristic  and  the  increases  in  the  cycle  length 
at  each  step,  and  therefore  the  total  number  of  iterations, will  be  the 
same  in  both  heuristics. 

It  should  also  be  noted  that  if  the  two  algorithms  have  the  same  values 
for  the  cycle  length,  then  the  values  for  the  best  unrestricted  n^'s 
found  using  (2.26)  at  each  iteration  will  also  be  identical.  By  examining 
a selected  group  of  problems  for  which  the  original  heuristic  performed 
worse  than  the  Goyal  procedure,  (Table  2.2)  it  was  found  that  if  in 
addition  to  examining  the  largest  lattice  point  neighboring  the  non-integer 
vector  of  n^'s  we  also  examined  the  closest  lattice  point,  then  the 
optimal  solution  remained  among  those  solutions  examined.  This  rule 


. Selected  Examples : Heuristic  Performs  Worse  than  Goyal 

Table  2.2 

r I 
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performed  better  than  picking  n*(T)  using  (2.15). 

One  can  also  mathematically  justify  the  use  of  the  closest  lattice 
point.  Let  us  consider  a two  retailer  problem  for  the  moment.  Given 
a particular  value  for  the  cycle  length  T,  the  cost  of  the  single 
cycle  policy  where  retailer  1 produces  n^  times  between  each  warehouse 
production  and  retailer  2 n^  times  can  be  written: 


fCn^njT)  = + 


Th0D0  nlKl  ThlDl  n2K2  Th2D2 

2 T 2n  T 2n2 


(2.32) 


the  gradient  of  f(n^,n2|T)  can  be  written 


/K  Th  D K Th  D 

"'■rtlT'jf-r'f-r 

V 2nl  2n2  . 


(2.33) 


The  unrestricted  (n^,n2)  pair  is  chosen  by  (2.26)  which  is  equivalent 
to  setting  7f(n^,n2|T)  = 0.  Let  us  approximate  the  cost,  given  T,  at 
the  lattice  point  (n^,n2)  by  expanding  the  cost  function  about  (n^,n2). 
That  is : 


f(n^,n2|T)  = f(n^,n^jT)  + (n*-n“,  n^-n^^ 7f(n^,n2|T) 


t 1,  i U i U._2£/  U Ui-w  l u £ U.t 
+ 2<ni~ni»  n2"n2)7  f(n1,n2|T)(n1-n1,  n2-n2)  . 


(2.34) 


By  definition,  7f(n^,n2|T)  = 0,  and  from  (2.33) 


72f(n^,n2|T)  = 


ThlDl 

(n^)3 


Th„D„ 
2 2 

, u. 3 

n2  . 


(2.35) 
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Therefore,  we  can  rewrite  (2.34)  as 


et  l | m \ c,  U Ui.v  ^ 1.  l U.2  ( 1D1  \ t 1.  I U 

f(n1,n2|T)  = f(nx,n  |T)  + j(n  -n  ) I — jj-jl  + ^n2_n2 

Vnp  J 


(2.36) 


and  clearly  the  minimum  functional  increase  in  (2.35)  over  the  value  of 
f(n^,n2|T)  will  take  place  at  the  closest  lattice  point  to  (n^n^). 
Therefore,  we  have  demonstrated: 


Lemma  2.2.  If  we  have  computed  the  best  unrestricted  n^(T),  i = 1,...,N, 

for  a given  value  of  T,  then  the  closest  lattice  point  to 

(n^(T) , . . . ,nN(T) ) will  have  the  smallest  increase  in  cost  compared  to 

the  cost  at  (n,  (T) , . . . ,n..(T) ) . □ 

1 N 

To  summarize,  the  revised  heuristic  will  examine  only  2 lattice  points 

neighboring  the  unrestricted  n^(T)'s  at  each  iteration:  the  closest 

c c 

lattice  point,  denoted  (n^(T) , . . . ,n^(T) ) and  the  largest  lattice  point, 
denoted  (n, (T) , . . . ,n„(T) ) . To  review,  let  us  rewrite  the  revised  heuristic 
as  follows : 


Step  0:  Set  the  active  node  equal  to  (1,1, ...,1).  Set  the  incumbent 
solution  equal  to  (1,1,..., 1)  with  cost  C*(l,l, . . . ,1) 
using  2.9.  Compute  T,  the  upper  bound  on  T*,  using  (2.28). 
Go  to  step  1. 

Step  1:  Compute  T = T*  (active  node).  If  T > T,  halt.  Otherwise 
go  to  step  2. 

Step  2:  Compute  the  best  unrestricted  n^(T)  (insuring  n^(T)  >_  1) 
i = 1,...,N  using  (2.26).  If  (n?(T) , . . . ,n^(T) ) , the 
closest  lattice  point,  was  not  previously  evaluated , then 
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compute  the  cost  and  compare  it  to  the  incumbent.  Go  to 
step  3. 

Step  3:  If  (n,  (T) , . . . ,n.,(T) ) , the  largest  lattice  point,  was  not 
1 N 

previously  evaluated  then  compute  its  cost,  compare  this 
cost  to  the  incumbent,  set  the  new  active  node  equal  to 

(n,  (T) , . . . ,n..(T) ) , and  go  to  step  1;  otherwise,  halt. 

IN 

2 . 10  Computational  Comparisons  for  Revised  Heuristic 

The  decision  rules  in  the  revised  heuristic  were  based  in  part 
upon  the  performance  of  the  original  heuristic  on  a particular  series 
of  test  problems.  To  asse?^  the  performance  of  the  revised  heuristic 
adequately  in  comparisor  to  the  Goyal  algorithm,  it  was  necessary  to 
generate  a new  series  of  175  test  problems.  The  problems  ranged  in 
si7.e  from  3 to  7 retailers,  the  set-up  costs  were  uniformly  distributed 
on  the  integers  from  1 to  100,  and  the  demand  rates  were  uniformly 
distributed  on  the  integers  from  1 to  10.  The  holding  costs  were 
uniformly  distributed  on  the  integers  from  1 to  10,  from  1 to  100  and 
from  1 to  1000  for  various  groups  of  problems.  The  results  for  these 
175  test  problems  are  summarized  in  Table  2.3.  In  these  175  test 
problems,  the  Goyal  algorithm  averaged  8.43  comparisons  to  the  incumbent 
per  problem  and  the  revised  heuristic  only  2.22  comparisons,  a 73.6% 
improvement.  In  addition,  even  for  these  relatively  small  test  problems, 
the  revised  heuristic  solved  the  175  problems  in  3.16  seconds,  the 
Goyal  algorithm  in  3.5  seconds.  For  these  small  problems,  the  average 
computation  time  is  still  very  short  for  both  procedures. 
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No.  of 
Retailers 

No.  of 
Problems 

Holding 

Cost  Range 

# Cor 
Goyal 

iparisons 

Heuristic 

# Optimal 

3 

25 

1-10 

109 

31 

25 

3 

25 

1-100 

121 

31 

25 

3 

25 

1-1000 

180 

45 

25 

5 

25 

1-10 

119 

56 

24 

5 

25 

1-100 

342 

79 

23 

5 

25 

1-1000 

300 

63 

24 

7 

25 

1-1000 

304 

83 

25 

Totals : 

175 

1-1000 

1475 

388 

171 

Performance  of  Revised  Heuristic  vs.  Goyal  Algorithm 

Table  2.3 

The  examination  of  only  2 lattice  points  has  eliminated  the  problem 
of  the  original  heuristic  examining  too  many  possible  solutions  in  a 
significant  proportion  of  the  problems.  In  only  1 example  out  of  175  does 
the  revised  heuristic  involve  more  comparisons  than  the  Goyal  procedure 
and  the  difference  is  only  a single  comparison.  (For  this  one  example 
the  heuristic  reevaluates  a previously  examined  solution  because  it 
only  checks  this  solution  for  repetition  against  those  generated  at  the 
preceding  iteration  rather  than  2 iterations  back.) 

Although  the  revised  heuristic  requires  less  computational  effort 
than  the  Goyal  approach,  it  no  longer  finds  the  optimal  solution  for 
every  test  problem.  In  four  out  of  175  cases  (Table  2.4)  or  2.3% 
the  heuristic  fails  to  find  the  optimal  solution.  Essentially,  in 
these  examples  the  heuristic  is  examining  the  wrong  lattice  point  at  a 
particular  iteration.  The  average  absolute  error,  for  these  four  problems 


i 


f 


46 


No.  of 
Comparisons 

To  Incumbent 

Cost  of 

Final  Incumbent 

Solution 

Absolute 

Error 

% 

Error 

Goyal 

Revised 
Heur . 

Goyal 

Heur. 

Goyal 

Heur. 

20 

6 

771.9117 

772.2543 

(1,1, 1,4, 6) 

(1,1, 1,3, 6) 

.3426 

.044% 

75 

9 

757.9408 

758.6303 

(1,1, 2, 6, 7) 

(1,1, 2, 5, 6) 

.6895 

.091% 

11 

6 

3378.9765 

3382.9260 

(.1,1, 1,1,4) 

(1,1, 1,1, 3) 

3.9495 

.12% 

12 

5 

369.4374 

379.0333 

(1,2, 2, 2, 3) 

(1,2, 2, 3, 3) 

.5959 

.16% 

Examples  Where  Revised  Heuristic  Fails  to  Find  Optimum 

Table  2.4 


in  error,  is  1.39  and  the  average  percentage  error  is  only  .1%,  certainly 
small  enough  to  be  acceptable  for  all  practical  purposes. 

The  improved  heuristic  was  also  tested  on  a series  of  larger  problems 
with  real  demand  rates  and  holding  costs  for  a large  corporation.  There 
were  20  retailers  in  these  problems  each  with  identical  set-up  and 
holding  costs.  The  demand  rates  for  retailers  2-20  were  integer 
multiples  of  the  demand  rate  at  retailer  1 ranging  from  2D  to  lOOD^. 
The  holding  cost  at  the  warehouse  was  1.06  times  the  holding  cost  at 
the  retailer  and  the  set-up  cost  at  the  warehouse  was  varied  to  test 
the  effectiveness  of  the  two  procedures.  The  results  for  the  two 
algorithms  are  summarized  in  Table  2.5.  Because  of  the  similarity  in 
cost  structure,  it  is  only  necessary  to  examine  selected  examples 
(Theorem  2.5). 

In  all  the  examples  tested,  the  heuristic  found  the  optimal  solution. 
More  importantly,  the  heuristic  seemed  to  remain  as  effective,  both  in 
terms  of  the  number  of  comparisons  and  in  time,  as  the  set-up  cost  at 


2K1 

7 

i 

5 

i 

.07 

.07 

10K 

7 

4 

.06 

.07 

100^ 

40 

8 

.13 

.08 

1000K 

193 

6 

. 

> .41 

.08 

Selected  Computational  Comparisons  for  Revised  Heuristic  vs. 

Goyal  on  20  Retailer  Problems 

Table  2.5 

the  warehouse  was  increased  relative  to  the  set-up  cost  at  the  retailers. 

The  Goyal  procedure,  on  the  other  hand,  became  much  worse  than  the  revised 

heuristic,  by  either  measure,  as  the  set-up  cost  increased.  As  the 

set-up  cost  at  the  warehouse  increases,  the  number  of  times  the  retailers 

produce  within  each  cycle  increases  as  well.  The  Goyal  procedure  only 

increases  one  n.  by  one  at  each  iteration  and  the  heuristic  can  make 
1 

larger  jumps  and  save  a great  deal  of  computational  effort. 

Graves  and  Schwarz  (1977)  have  presented  a branch  and  bound  algorithm 
to  find  the  optimal  single  cycle  policy.  Although  no  direct  comparisons 
have  been  made,  the  improved  heuristic  seems  superior  for  larger  numbers 
of  retailers.  There  have  been  no  indications  that  the  improved  heuristic 
fails  to  find  the  optimal  solution  as  frequently  for  problems  with  a 
larger  number  of  retailers  as  was  the  case  for  the  branch  and  bound 
algorithm  of  Graves  and  Schwarz.  Although  computation  time  seems 
comparable  for  the  two  approaches  on  small  test  problems,  the  enumeration 
tree  in  the  branch  and  bound  algorithm  has  a level  for  each  retailer  and 
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for  extremely  large  problems,  the  storage  requirements  may  be  unwieldly. 
Also,  at  every  level  of  the  enumeration  tree  there  is  a branch  for  each 
integer  value  of  n^.  Thus,  for  those  problems  where  the  n/s  get 
large  and  the  Goyal  algorithm  performs  much  worse  than  the  revised 
heuristic,  the  enumeration  tree  in  the  branch  and  bound  algorithm  also 
will  be  wide  at  each  level.  Although  no  definitive  statements  can  be 
made  without  further  direct  comparisons,  the  iterative  heuristic  seems 
at  least  as  effective,  if  not  more  so  for  large  problems,  than  the 
existing  solution  techniques. 

2 . 11  Single  Cycle  Policies  and  Optimality 

With  the  computational  techniques  discussed  in  sections  (2. 6)- (2. 10) , 
it  is  possible  to  determine  the  optimal  single  cycle  policy.  The  major 
remaining  question  is  whether  a single  cycle  policy  is  the  overall 
optimal  policy.  The  answer  to  this  question  is  yes,  in  the  cases  of 
one  retailer  or  N identical  retailers  from  theorems  2.1  and  2.2,  but 
we  are  interested  in  examining  the  correlation  between  single  cycle 
policies  and  overall  optimal  policies  in  more  general  cases. 

Let  us  begin  by  recalling  the  example  used  to  demonstrate  the 
heuristic  algorithm: 


Kx  = 99.9,  h1  = 99,  D,  = 1,  (2.31) 

K2  = 99.9,  h2  = 199,  and  P2  = 1. 


For  this  example,  the  optimal  single  cycle  policy  has  (n^.nj)  = (2,3) 
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with  an  average  cost  per  unit  time  of  343.13125.  Let  us  split  the  problem 
into  two  independent  one  warehouse,  one  retailer  problems,  and  then  deter- 
mine the  best  "separate  retailing"  policy  as  discussed  in  Section  2.5. 

F irst , we  compute 

KQh1/K1h0  = 9.9/99.9 

and  KQh2/K2h0  = 19.9/99.9 

and  in  both  cases,  the  smallest  integer  n such  that  n(n+l) 

than  these  quantities  is  one  and  thus  best  separate  retailing 
SR  SR 

(n1  ,n2  ) = (1,1)  and  the  average  cost  per  unit  time  of  this 
from  (2.18)  is : 

[2(K0+l-K1)(h0D1+h1D1/l)]1/2  + [2(K0+l-K2)(h0D2+h2D2/l)]1/2  = 341.42. 


is  greater 
policy  has 
policy 


We  have  therefore  demonstrated. 

Theorem  2.6.  The  overall  optimal  policy  may  not  be  a single  cycle  policy 
and 


Corollary  2.7.  The  best  separate  retailing  policy  may  be  less  costly  than 
the  best  single  cycle  policy. 

In  the  best  separate  retailing  policy  described  above,  the  retailers 
have  optimal  lot  sizes  of 


/2  and 


(2.37) 


1 
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where  and  T0  can  be  derived  from  (2.13)  for  one  retailer.  This 

separate  retailing  policy  is  stationary  at  the  retailers  but  not  at  the 
warehouse.  The  warehouse  produces  lots  of  size  ^2  when  it  is  producing 
for  retailer  1 and  lots  of  size  1 when  it  is  producing  for  retailer  2. 

Consider  now  the  stationary  policy  based  upon  the  separate  retailing 
plan  described  above.  The  retailers  again  have  lot  sizes  of  J2  and  1; 
but , the  warehouse  produces  a lot  size  of  1 + ^2  whenever  its  on-hand 
inventory  is  insufficient  to  meet  a retailer's  demand  requirements.  To 
evaluate  the  average  cost  per  unit  time  for  this  stationary  policy , it 
should  be  noted  that,  on  the  average  the  warehouse  sets-up  once  every 
(1  + /2)/2  time  units  and  the  average  on-hand  inventory  per  unit  time 
at  the  warehouse  is  (1  + /2)/2.  Thus  the  average  cost  per  unit  time 
for  this  policy  is: 


+ 200(j) 


warehouse  costs  retailer  retailer  holding 

set-up  costs  costs 


342.54. 


Thus  this  stationary  policy,  although  more  costly  than  the  best  separate 
retailing,  is  less  costly  than  the  best  single  cycle  policy.  □ 

The  stationary  policy  described  above  cannot  be  a single  cycle  policy 
since  there  are  no  renewal  points  after  time  zero.  (Due  to  the 
irrationality  of  /T,  retailers  1 and  2 never  produce  simultaneously 
after  time  0).  This  example  therefore  proves. 
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Theorem  2.8.  Among  the  class  of  stationary  policies,  a single  cycle 
policy  need  not  be  optimal.  □ 

Theorem  2.8  contradicts  Theorem  1 of  Graves  and  Schwarz  ( 1977 ) which 
states,  that  the  optimal  stationary  policy  for  a one-warehouse,  N 
retailer  system  is  indeed  a single  cycle  policy.  Their  proof  of 
Theorem  1 is  in  error  because  it  relies  upon  the  application  of  proper- 
ties PI  and  P4  to  the  optimal  stationary  policy.  These  properties  which 
hold  for  the  overall  optimal  policies  need  not  be  valid  for  optimal 
stationary  policies. 

Even  if  we  had  not  previously  computed  the  cost  of  the  best  separate 
retailing  policy,  it  is  obvious  that  the  stationary  policy  constructed 
above  cannot  be  the  overall  optimal  solution.  The  irrationality  of 
✓2,  in  addition  to  precluding  the  existence  of  renewal  points  as 
discussed  above,  also  forces  the  warehouse  to  produce  when  its  on-hand 
inventory  is  greater  than  zero.  This  violates  property  PI  which  holds 
for  all  overall  optimal  policies.  Graves  and  Schwarz  (1978)  have 
suggested  that  if  we  restrict  our  attention  to  the  class  of  stationary 
policies  which  satisfy  P1-P5,  then  among  this  class  of  policies,  a 
single  cycle  policy  is  optimal.  This,  however,  is  also  not  the  case  and 
we  can  prove 

Theorem  2.9.  Among  the  class  of  stationary  policies  satisfying  P1-P5, 
a single  cycle  policy  need  not  be  optimal. 

Proof  of  Theorem  2.9. 


The  proof  will  be  by  construction  of  an  example  for  which  there 
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exists  a stationary  policy  satisfying  P1-P5  with  lower  cost  than  the  best 
single  cycle  policy. 

Consider  the  following  one-warehouse,  two  retailer  example: 


K.  = .1,  h = 1,  D = 1 + /2, 

0 ’ 0 ’ 0 ’ 

Kx  = 34.9,  h = 70/2  - 1,  D = /2, 


K2  = 99.9,  h2  = 199,  D2  = 1. 


(2.38) 


For  this  example,  one  can  determine  that  the  best  single  cycle  policy 

has  (n^,n2)  = (3,2)  with  average  cost  per  unit  time  equal  to  300.38022. 

SR  SR 

The  best  separate  retailing  policy  for  this  problem  has  (n^  ,n  ) = 
(1,1)  and  an  average  cost  per  unit  time  of  298.99494.  The  optimal  lot 
sizes  for  this  separate  retailing  policy  are: 


(2.39) 


Therefore  retailer  1 produces  a lot  of  size  1 every  1//2  time 
units,  retailer  2 produces  a lot  of  size  1 every  1 time  unit  and 
the  warehouse  produces  a lot  of  size  1 each  time  either  of  the 
retailers  produce,  except  at  time  zero  where  the  warehouse  produces  2. 
Let  us  examine  the  stationary  policy  where  the  retailers  follow  the 
same  policy  as  in  separate  retailing  and  the  warehouse  produces  2 
each  time  it  produces.  The  average  time  between  productions  at  the 
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warehouse  is  2/DQ  = 2/(1  + ^2)  and  the  average  inventory  held  physically 
at  the  warehouse  is  1/2.  The  average  cost  per  unit  time  for  this 
stationary  policy  is: 


C 


s 


1(i)  + (34-9)(,/2)  + 99-9  + 100  + 35/2 

warehouse  costs  retailer  set-up  retailer 

costs  holding  costs 


= 299.37423. 


Due  to  the  irrationality  of  /2,  there  is  no  point  at  which  retailer  1 
and  retailer  2 produce  simultaneously  after  time  0.  Thus  there  are  no 
points  after  time  zero  where  the  system  inventory  is  zero,  hence  no 
renewals.  Due  to  the  fact  that  both  retailers  have  identical  lot  sizes, 
the  warehouse  will  only  produce  when  it  has  zero  inventory  on-hand.  Thus, 
property  PI  is  satisfied  for  this  stationary  policy.  Properties  P2-P5 
are  obviously  also  satisfied  and  therefore  this  policy  is  stationary, 
satisfies  P1-P5,  and  has  lower  cost  than  the  best  single  cycle  policy 
thereby  establishing  Theorem  2.9.  □ 

As  mentioned  earlier,  in  the  proof  of  their  Theorem  1,  Graves  and 
Schwarz  (1977)  applied  properties  PI  and  P4  to  the  optimal  stationary 
policy  which  may  not  be  valid.  However,  if  we  restrict  our  attention 
to  the  class  of  stationary  policies  satisfying  properties  PI  through 
P5,  then  the  proof  of  Theorem  1 would  seem  to  demonstrate  that  among  this 
class  of  policies,  a single  cycle  policy  is  optimal.  Theorem  2.9  disproves 
this  conjecture  and  therefore  there  must  be  an  additional  flaw  in  the 
arguments  presented  by  Graves  and  Schwarz.  The  flaw  lies  in  the 
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assertion  that  stationarity  and  property  P4  together  imply  that  there 
exists  an  integer  n^  such  that  every  n^th  production  run  at  the 
warehouse  is  a simultaneous  production  point  for  retailer  k.  Although 
property  P4  implies  that  the  time  interval  between  simultaneous  produc- 
tion points  is  finite,  it  is  not  necessarily  the  case  that  this  time 
interval  is  also  constant.  In  the  example  used  to  prove  Theorem  2.9, 
the  irrationality  of  and  the  lack  of  renewal  points  causes  the 

time  interval  between  simultaneous  production  points  to  differ  over 
time.  Thus,  there  is  no  single  integer  n^  such  that  every  n^th 
production  point  at  the  warehouse  is  also  a production  point  for 
retailer  k. 

To  summarize,  we  have  demonstrated  in  Theorems  2.6,  2.8  and  2.9 
that  a single  cycle  policy  need  not  be  optimal  for  increasingly 
more  restrictive  classes  of  policies.  A common  property  of  all  the 
examples  presented  in  this  section  is  the  lack  of  renewal  points 
in  those  policies  that  are  less  costly  than  the  best  single  cycle 
policy.  Further  investigation  of  single  cycle  policies  as  a member 
of  the  class  of  renewal  policies  is  warranted  and  this  is  the  subject 
of  the  next  section. 

2.12  Single  Cycle  Policies  and  Renewal  Policies 

A single  cycle  policy  is  a special  case  of  a renewal  policy. 

In  this  section,  we  will  derive  results  concerning  the  holding  cost  in 
multiple-cycling  policies  and  then  demonstrate  that,  in  certain 
instances,  single  cycle  policies  are  indeed  optimal. 

We  begin  by  noting  that  the  example  used  to  prove  Theorem  2.9  has  a 
stationary,  non-renewal  policy  that  is  less  costly  than  the  best  single 
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cycle  policy.  The  following  theorem  relates  optimal,  stationary  and 
renewal  policies. 


Theorem  2.10.  If  the  overall  optimal  policy  is  stationary,  then  it  is  a 
renewal  policy. 


Proof  of  Theorem  2.10. 

We  begin  by  assuming  that  the  optimal  policy  is  stationary  and  has 
no  renewals.  In  this  case,  we  will  construct  a non-stat ionary  policy 
which  has  lower  average  cost  per  unit  time  than  the  hypothesized  optimal 
policy. 

By  stationarity  and  property  PI,  retailer  i has  a production  point 
every  T^  = Q^/D^  time  units.  By  property  P2,  the  warehouse  produces 
only  when  at  least  one  of  the  retailers  produces  and  thus , without  loss 
of  generality,  we  may  assume  that  the  first  warehouse  production  after 
time  zero  is  a simultaneous  production  point  for  retailer  1;  that  is, 
the  first  warehouse  production  point  after  time  zero  is  at  time  k T , 

rkiTii 

integer.  During  k^T^,  — - — lots  of  size  , including  the 

lot  at  time  zero,  are  required  for  retailer  j,  where  fxl  is  the 


- r^i 

3 f . | • 


smallest  integer  greater  than  or  equal  to  x.  Let  k 
klTl 

-= — cannot  be  an  integer  for  all  j i 1,  or  k T would  be  a 

j 

renewal  point,  contradicting  the  hypothesis. 

Since  property  PI  must  be  satisfied  for  the  optimal  solution,  the 
warehouse  must  have  zero  inventory  on-hand  at  k.^.  Therefore, 


N fk  T“|  N 

*0  ■ -A  + ■ kA  Ma  Vi- 


(2.41) 


56 


Since  we  have  assumed  that  the  first  warehouse  production  occurs  at  k^T^, 
it  is  clear  that 

(k.-l)T.  < k,T,  < k.T.  for  j t 1.  (2.42) 

: : ii-:: 

In  fact,  this  inequality  must  be  a strict  inequality  for  at  least  one 
retailer  or  the  system  would  have  a renewal  at  time  k^T^. 

It  is  also  clear  from  (2.41)  and  stationarity  that  the  length  of 
time  between  successive  warehouse  production  points  is  greater  than 
or  equal  to  min{k^T  ,k  T2 , . . . ,^1^} , which  is  equal  to  k^T^  by  (2.42). 
Thus,  all  warehouse  production  intervals  have  length  at  least  k^T  . 

Now  consider  those  intervals  where  the  on-hand  inventory  is  zero 
at  the  warehouse.  Due  to  the  lack  of  renewals,  intervals  with  lengths 
equal  to  infinitely  many  values  between  0 and  min(T^,T  ,. . . ,TN)  may 
be  found.  Let 

6 = k T -min(Q  /(Q0+Qi)).  (2.43) 

i 

Look  at  an  interval,  (p1»?2),  with  zero  on-hand  inventory  at  the 
warehouse  beginning  at  time  p^  and  lasting  until  time  p2  and  length 
6',  6'  < 6.  Let  pQ  be  the  time  of  the  warehouse  production  point 

immediately  preceding  p^.  It  should  be  noted  that  pQ  and  p2  are 
warehouse  production  points,  while  p^  is  a production  point  for  a 
retailer,  or  some  group  of  retailers. 

In  the  present  policy,  Qq  units  are  produced  at  the  warehouse 
at  times  pQ  and  p2,  and  zero  units  are  produced  at  the  warehouse 
at  time  p . At  least  one  retailer  produces  at  time  p^,  and  we  will 
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denote  the  total  produced,  at  this  time,  over  all  retailers  by  Q.  Clearly, 


Q ^ Q_.  for  some  j . (2.44) 

Consider  a policy  which  is  identical  to  the  original  optimal  policy 
at  all  the  retailers  but  is  revised  at  the  warehouse  as  follows: 


(1) 

Produce 

o 

O' 

- Q 

at 

time 

(2) 

Produce 

Qo 

+ Q 

at 

time 

(3) 

Produce 

0 

at  time 

P2 " 

The  only  cost  differences  from  the  original  production  plan  are  over 
the  interval  (pQ,P2).  The  set-up  costs  remain  unchanged;  but,  in  the 
revised  policy,  Q fewer  units  are  held  over  the  time  interval  (Pq’P]^ 
and  Qq  additional  units  are  held  over  the  time  interval  (p^»P2)  as 
compared  to  the  original  plan.  Therefore,  the  net  change  in  inventory 
of  this  revision  equals 


Wpi}  ■ Q(prpo): 


Let  us  assume  that 

6'  = (P2_P1)  - (Q/Qo^pi"po)i  (2.45) 

that  is,  the  revised  plan  carries  at  least  as  much  inventory  at  the  ware- 
house as  the  original  plan. 

As  discussed  earlier,  all  warehouse  production  intervals  have  length 
at  least  k^T^  and  therefore 
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p2  - p0  - W 


(2.46) 


But 


P2  " P0  = (P2_P1)  + (P1_P0)  = 6'  + (pi-p0>«  From  (2.46), 


P1  - P0  iklTl  - s' 


and,  therefore,  (2.45)  implies 


6'  1 (Q/Q0)(k1T1-6'). 


(2.47) 


Since 


Q > Qj 


for  some  j , 


<$'  (Q_. /Qq ) (k^T^-6 ' ) for  some  j. 


Rearranging  terms,  we  see  that 


/Qo+Qi\  Qi 

6'  ^ — q— M >_  ^4-  (k1T1 ) for  some  j, 


or 


"(^)klTl"klTl’1n(vQi, 


= 6 by  definition. 


This  is  clearly  impossible  by  the  choice  of  <5 ' , which  is  strictly  smaller 
than  5.  Therefore,  this  revised  plan  has  fewer  units  of  inventory  carried, 
and  hence  lower  cost,  than  the  original  over  the  time  interval  (P0,P2^* 

By  repeating  this  construction  for  all  intervals  with  zero  on-hand 
inventory  at  the  warehouse  and  length  less  than  6 , we  have  demonstrated 
that  there  exists  a policy  with  lower  average  cost  per  unit  time  than  the 
hypothesized  optimal  policy.  This  establishes  Theorem  2.10.  □ 
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Theorem  2.10  indicates  that  further  examination  of  stationary  renewal 
policies  is  in  order.  In  stationary  renewal  policies,  a renewal  point 
will  occur  at  every  n^th  warehouse  production  point,  where  n^  _>  1 
and  integer.  Clearly,  if  nQ  = 1,  we  have  the  usual  single  cycle  policy. 
If  nQ  > 1,  we  will  call  this  policy  a multiple  cycle  policy. 

We  will  demonstrate  that  for  a particular  problem,  the  best  single 
cycle  policy  may  have  higher  average  cost  per  unit  time  than  a multiple 
cycle  policy.  This  is  not  as  obvious  a fact  as  it  may  seem.  Multiple 
cycle  policies  do  constitute  a wider  range  of  policies  than  single 
cycle  policies,  since  nQ  is  not  constrained  to  be  1.  But  the  warehouse 
echelon  inventory  in  a multiple  cycle  policy  may  not  be  easily  computed. 

To  illustrate  some  of  the  properties  of  a multiple  cycle  policy,  let 
us  consider  the  multiple  cycle  policy  with  (n^n^n^  equal  to  (2,2,3) 
and  T equal  to  the  time  between  renewals.  We  begin  by  developing  an 
expression  to  represent  the  average  cost  per  unit  time  of  this  policy. 

Retailer  1 produces  lots  of  size  at  times  0 and  T/2. 

Retailer  2 produces  lots  of  size  at  times  0,  T/3  and  2T/3. 

In  this  multiple  cycle  policy,  as  in  a single  cycle  policy,  retailer  1 
holds  an  average  of  Q^/2  units  per  unit  time,  and  retailer  2 holds 
an  average  of  Q^/2  units  per  unit  time. 

By  stationarity , the  warehouse  produces  two  lots  of  size  Qq,  where 


If  the  warehouse  was  to  produce  at  times  0 and  T/2,  then  the  average 
amount  of  inventory  held  per  unit  time  would  be  Qq/2  as  in  a single 
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cycle  policy.  Unfortunately,  this  plan  is  infeasible. 

Consider  Figure  2.4,  which  illustrates  the  amount  of  echelon  inventory 
in  the  multiple  cycle  policy  (2,2,3).  In  addition,  the  on-hand  inventory 
at  the  warehouse  is  represented  by  the  shaded  area.  At  time  zero,  Qq 
is  produced  at  the  warehouse  and  lots  of  and  are  sent  immediately 

to  retailers  1 and  2,  respectively.  From  (2.48),  this  leaves  the 
warehouse  with  an  on-hand  inventory  of  Q2/2.  Ar  time  T/3,  retailer  2 
requires  another  lot  of  size  Q2  from  the  warehouse.  But  the  on-hand 
inventory  is  insufficient  to  meet  this  need.  To  maintain  feasibility, 
the  warehouse  must  produce  a lot  of  size  at  time  T/3. 

Recall  that  if  the  warehouse  produces  at  times  0 and  T/2,  then 
the  average  echelon  inventory  is  Qq/2.  In  the  actual  plan,  the  ware- 
house must  produce  at  T/3  rather  than  T/2  and  thus  the  total  echelon 
inventory  for  the  cycle  T is  increased  by  Qq(T/2  - T/3)  = QqT/6. 

The  increase  in  the  average  echelon  inventory  per  unit  time  is  Qq/6. 

Thus  the  average  warehouse  echelon  inventory  per  unit  time  is 

V?  * ?>  * I V <2-u9) 


For  each  facility  i,  including  the  warehouse,  n^K^  is  the  cost 
of  set-ups  in  each  cycle  of  length  T.  Thus,  the  average  set-up  cost 
for  facility  i,  per  unit  time  is  n^K^/T. 

Given  (2.49),  we  can  now  express  the  average  cost  per  unit  time 
of  the  multiple  cycling  policy  (2,2,3)  as 


2K„ 


+ 2K  + 3K 

— r ~ + Vi  V + 


hi(iV  + h2(iQ2)- 


(2.50) 
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FIGURE  2.4  MULTIPLE  CYCLE  POLICY  WITH  nQ*2,  2,  n2»  3 
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By  stationarity , 
^n0,ni’n2^  ~ (2,2,3) 

2Kq  + 


= TD^/n^  for  all  i.  Then  for  this  example  with 
we  can  rewrite  (2.50)  as 


2KX  + 3K2 

T 


+ T 


(Vo  Vi 

\ 3 + 4 


(2.51) 


Therefore,  given  the  choice  of  (nQ,n^,n  ) = (2,2,3),  it  is  clear  that  the 
optimal  cycle  length  T*(2,2,3)  is  given  by 


T*(2,2,3) 


2KQ  + 2KX  + 3K3 

h0°0  , hlDl  , h2°2 
^ 3 + 4 + 6 - 


1/2 


(2.52) 


and  the  average  cost  per  unit  time  obtained  when  T*(2,2,3)  is  substituted 
into  (2.51)  is  given  by 


/h  D h D h D \ 

C*( 2,2,3)  = 2 (2KQ  + 2KX  + 3^)  + ~~  + 


1/2 


Consider  the  following  one  warehouse,  two  retailer  system  with  costs 
and  demand  rates  as  follows : 


Ko  = x*  h0  = D0  = 2’ 


Kx  = 17,  h = 10,  D = 1, 


K2  = !7,  h2  = 22,  and  D2  = 1. 


(2.54) 


We  can  compute  the  cost  of  the  multiple  cycle  policy  (2,2,3)  for  this 
example  using  (2.53)  and  we  find 


C*(2 ,2 ,3 ) = 48.7647. 
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However,  the  best  single  cycle  policy  (1,1)  has  cost 

, 

C*(l,l)  = 48.7852 

and,  therefore,  we  have  demonstrated 

Theorem  2.11.  Among  the  class  of  stationary  renewal  policies,  a single 
cycle  policy  may  not  be  optimal.  □ 

It  is  clear,  however,  that  the  multiple  cycling  policy  discussed  ■ 

above  cannot  be  the  overall  optimal  solution  either.  Note  that  property 

PI  is  violated  at  time  T/3.  That  is,  the  warehouse  produces  at  time  ■ 

T/3  with  an  on-hand  inventory  of  Q2/2.  Clearly,  if  the  warehouse 

produces  Q^/2  fewer  units  at  time  zero  and  Q2/2  more  units  at  time 

T/3,  feasibility  will  be  maintained  and  the  average  inventory  per  unit 

time  at  the  warehouse  will  be  lower  than  in  the  multiple  cycling  policy 

discussed  above.  Thus,  there  is  an  easily  constructed  non- stationary 

policy  with  lower  average  cost  per  unit  time  than  the  multiple  cycling 

policy  (2,2,3). 

Furthermore,  it  is  possible  to  perform  the  same  type  of  shifting 
of  production  for  any  general  multiple  cycle  policy,  thus  demonstrating 
that  a multiple  cycle  policy  cannot  be  the  overall  optimal  policy.  We 
know  by  Theorem  2.10,  that  if  the  overall  optimal  policy  is  stationary 
then  it  is  a renewal  policy.  Therefore,  we  have  shown 


Theorem  2.12.  If  the  overall  optimal  policy  is  stationary,  then  it  is 
a single  cycle  policy.  □ 
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We  would  also  like  to  compare  the  quality  of  single  cycle  policies 
to  multiple  cycle  policies  in  more  general  cases.  In  order  to  do  this, 
it  is  necessary  to  derive  a general  expression  for  the  echelon  holding 
costs  at  the  warehouse,  as  we  did  earlier  for  the  multiple  cycle  policy 
(2,2,3).  Unfortunately,  it  is  impossible  to  derive  an  exact  expression 
based  solely  on  the  policy  (n  ,n  ,n2 » • • . »n^).  This  can  be  seen  by  the 
following  result. 

Theorem  2.13.  In  a multiple  cycle  policy  (n  ,n^, . . . ,n  ) , the  echelon 
holding  cost  at  the  warehouse  is  dependent  upon  the  lot  sizes 

^1*^2’"  ' ’V' 

Proof  of  Theorem  2.13. 

Consider  the  multiple  cycle  policy  with  (n^n^n^)  = (2,3,5).  Assume 

3 5 

that  the  cycle  length  is  T.  By  stationarity , Qq  = J ^i  + J ^2‘  At  t^me 
T ii 

the  warehouse  has  — + — Q2  units  on  hand  and  units  are 

required  for  retailer  1.  If  < Q^,  then  the  warehouse  must  produce 
T 

Qq  units  at  time  — . If  Q >_  Q , then  the  warehouse  need  not  produce 

T 2T 

at  time  — but  it  will  have  to  produce  at  — . Thus  if  Q2  < Q , the 

average  echelon  inventory  at  the  warehouse  is  Q_(i-  + ^-)  and  if  Q < Q0 , 

the  average  echelon  inventory  at  the  warehouse  is  QQ(y  + yy)-  This 

demonstrates  Theorem  2.13.  □ 

It  is  clear  from  the  example  above,  that  for  nQ  = 2,  the  echelon 
inventory  at  the  warehouse  is  determined  by  which  retailer  has  a simul- 
taneous production  point  with  the  warehouse  after  time  zero.  Similarly, 
if  nQ  > 2 , then  the  echelon  inventory  at  the  warehouse  is  determined  by 
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the  sequence  of  retailers  having  simultaneous  production  points  with  the 
warehouse.  The  problem  of  determining  the  warehouse  echelon  inventory 
for  a general  sequence  of  retailers  appears  impossible  but  we  can 
analyze  a special  case. 

Theorem  2.14.  For  the  multiple  cycle  policy  (nQ ,n  , . . . ,n^) , if  each 
warehouse  production  point  is  a simultaneous  production  point  for  the 
same  retailer  i , then  the  average  echelon  inventory  at  the  warehouse , 

n^  is  not  an  integer  multiple  of  nQ  (2.55) 

n^  is  an  integer  multiple  of  nQ  greater 
than  n . (2.56) 


I , can  be  expressed  by 


1 no-^ 

I,VTt-H7>  if 


and  is  bounded  below  by 


n -1 

1 * Vi  * -§r>  if 

l 


Proof  of  Theorem  2.14: 

If  n^  is  not  an  integer  multiple  of  n^,  then  we  can  write 

ni  = kn0  + £ w^ere  k is  an  integer  and  e is  an  integer  strictly 

between  0 and  nQ.  There  are  nQ  production  intervals  to  be  considered. 

R 

The  length  of  each  of  these  production  intervals  is  either  — T or 

k+1  ^ 

T,  where  T is  the  length  of  time  between  the  overall  renewals. 

If  there  are  b intervals  of  length  — T and  c intervals  of 

n . 

. k+1  1 

length  , then 

i 


k k+1 

b(— )T  + c(— )T  = T and  b+c  = n. 
n . n . 

l l 
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It  is  clear  that  b = nQ-e  and  c = e.  Furthermore,  the  ending  echelon 
inventory  for  each  of  the  warehouse  production  intervals  takes  on  the 

1.. ^ r\  i (C“l)  « -»_1 1 


values  — Q. , — Q. , . . . , £ 1 for  the  e production  intervals 

n.  0’  n.  v0’  ’ n.  v0  ^ 

11  1 , . (n  -1) 

with  length  — — T and  the  values  — Q , — — — Q , , Q for 

e n.  n.  0 n.  0 n.  0 

l l i l 

the  (nQ-e)  production  intervals  with  length  — T.  This  is  clear 
because  there  are  no  renewals  and  thus  each  value  between  — Q.  and 

X T v n*  1 

/n  -1\  l 

( I Q.  must  be  taken  on.  Also,  for  each  production  interval  that 

vni/  1 . 

lasts  for  — T time  units  the  echelon  inventory  is  ( — )Q.  units  higher 
n,  n . 0 

^ 1 

at  the  end  than  at  the  start  and  therefore  each  — T type  interval  must 

e n;L 

have  ending  inventory  of  at  least  — Q^. 

i 

Also,  since  it  requires  T/nQ  time  units  to  use  a lot  of  size  Q^, 
the  amount  of  inventory  used  in  an  interval  of  length  k/n^  is 
Q0(knQ/n^).  Similarly,  the  amount  of  inventory  used  in  an  interval  of 
length  (k+l)/n^  is  Q0((k+l)nQ/n^).  Therefore  we  can  write  the  total 
warehouse  echelon  inventory  for  a cycle  as 


£ (Average  height  of  production  interval) 

All  production 
intervals 

•(Length  of  production  interval) 


r , Inventory  used  in  production  interval  _ , . . . x 

l ( *■ + Ending  inventory) 

All  production 

intervals  with 


length  Tk/n^ 


•(Tk/n.) 

i 


c , Inventory  used  in  production  interval  _ , 

+ l ( *■ ^ + Ending  inventory) 

All  production 

intervals  with 


length  T(k+l)/n^ 


.(T(k+l)/n.) 


0 \ n ^ J 2n^  V n n^  ni  / ni 

* ■ ^ <V^>  * vY»  — V^<) 


(n0-e)( 


■ <v‘>f  »0{^i}  ♦ 


since  n.  = kn„+e  and 

l 0 


■{^) 


TQ, 


Therefore,  the  average  echelon  inventory  held  at  the  warehouse  per  unit 


time  can  be  expressed  as 
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If,  on  the  other  hand,  n^  = knQ  where  k is  an  integer  and 

greater  than  1,  then  the  first  warehouse  production  after  time  zero 

k-1 

occurs  no  later  than  at  time  T( ).  This  warehouse  production  cannot 

T ^ 

occur  at  time  — or  the  system  would  have  a renewal  at  that  point  in 

time.  The  ending  echelon  inventory  for  this  warehouse  production  interval 

k-1  /t  1 

is  at  least  as  large  as  Q_ ( 1 T/ — ) = Q„(t).  In  fact,  because  n. 

0 ni  / n0  0 k 1 

is  a multiple  of  nQ,  the  ending  echelon  inventory  cannot  be  any  lower 

than  Qg(j“)  until  the  final  warehouse  production  interval  in  the  cycle. 

k-1  k k k k+1 

If  the  warehouse  production  intervals  have  lengths  , — , — , ..., — , , 


-j 


then  the  ending  echelon  inventory  for  all  the  warehouse  production 
intervals,  except  the  last,  is  exactly  QgC^-).  So,  clearly,  the  average 
echelon  inventory  for  this  plan  is  a lower  bound  for  the  average  echelon 
inventory  in  the  actual  production  plan. 

If  we  compute  the  average  inventory,  I,  for  this  plan  we  find 


k 2 k 2 


Ending  echelon  Average  inventory  used  for 

inventory  each  warehouse  production  interval 


There  do  exist  some  special  cases  of  multiple  cycle  policies  where 
all  warehouse  production  points  must  be  simultaneous  production  points  for 
a particular  retailer. 
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Corollary  2.15.  The  formula  (2.55)  may  be  used  to  exactly  express  the 
average  echelon  inventory  at  the  warehouse  for  the  multiple  cycle  policy 

(no’ni’‘ " ‘ ,nu ^ ^ 

(a)  n^  ^ nQ,  for  some  i,  and  n^  = nQ  for  all  j i i or 

(b)  n.  = n^  for  all  i,j  not  equal  to  zero. 


Proof  of  Corollary  2.15.  * 

It  is  only  necessary  to  show  that  Theorem  2.14,  for  n^  not  a multiple 
of  n^,  may  be  applied.  In  both  cases  (a)  and  (b),  all  warehouse  produc- 

\ 

tions  are  triggered  by  a retailer  i production.  In  either  case,  if 

n^  = knQ,  then  all  the  n^  are  divisible  by  nQ  and  the  multiple  cycle 

°1  n2  nN  ' 

policy  is  equivalent  to  the  single  cycle  policy  (1,  — , — ,..., — ).  □ 

n„  n„ 


Suppose  we  evaluate  the  formulas  (2.55)  and  (2.56)  for  each  of  the 

retailers  in  the  multiple  cycle  policy  (n. ,n. , . . . ,n„) . Denote  the 

0 1 N 

minimum  of  these  quantities  by  L(n  ,n  , . . . ,n.  ) . For  some  special  cases 

U X N 

of  multiple  cycle  policies  L(n. ,n. ,. . . ,n„)  is  a lower  bound  on  the 

U 1 N 

average  echelon  inventory  held  at  the  warehouse. 


Corollary  2.16.  L(nQ ,n^, . . . ,n^)  is  a lower  bound  on  the  average  echelon 
inventory  at  the  warehouse  if 

(a)  nQ  = 2 or 

(b)  n^^  > n0n^  for  j i i. 

Proof  of  Corollary  2.16. 

(a)  If  nQ  = 2,  both  warehouse  production  points  are  simultaneous 
production  points  for  a particular  retailer.  Without  additional  information 
concerning  the  lot  sizes,  it  is  impossible  to  determine  which  retailer 
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should  be  used  to  evaluate  (2.55)  or  (2.56).  By  evaluating  (2.55)  or 
(2.56)  for  all  the  retailers  and  choosing  the  minimum,  we  must  have  a 


lower  bound  on  the  actual  average  echelon  inventory. 

(b)  If  n.  > n„n.  then  — < — — and  clearly  L(n„ ,n_ , . . . ,n„) 
l On  n.  n.n.  J 01  N 

J l 0 ] 

is  found  by  evaluating  (2.55)  or  (2.56)  for  n..  Suppose  the  (£+l)st 

1 k. 

warehouse  production  is  triggered  by  retailer  j at  time  — T where 

k. 

is  an  integer.  Then  the  total  echelon  inventory  is  CL( — T - — T) 

0 n0  ’ nj 

higher  for  this  plan  than  if  the  warehouse  produced  exactly  at  time 

k . 

all  integers  , T( — ^-)  > 


— T.  Since  l.  n. , k.  and  n.  are 
n0  0 3 3 

11 

Since  — < , by  assumption,  there  is  a retailer  i production 

ni  Vj 


n.  n . 

0 3 


Vj 


-i  £T 

between  — ^ T and  — and  the  total  echelon  inventory  would  be  lower  if 

nj  no 

it  were  possible  to  delay  the  (£+l)st  warehouse  production  point  and 
make  it  a simultaneous  production  point  for  retailer  i. 

This  same  argument  will  hold  for  all  warehouse  production  points  and 
therefore  the  possibly  infeasible  plan  where  all  warehouse  production 
points  are  simultaneous  production  points  for  retailer  i has  average 
echelon  inventory  no  larger  than  the  actual  plan.  Since  n^  minimizes 
L(nQ  ,n^, . . . ,n^) , this  is  clearly  a lower  bound  on  the  average  echelon 
inventory.  □ 


In  general,  it  is  possible  to  compute  a lower  bound  on  the  average 
echelon  inventory  for  the  multiple  cycle  policy  (nQ ,n  , . . . ,nN)  that 

is  better  than  the  obvious  lower  bcund  of  Qn/2.  Simply  determine  the 

T 2T  (V1)T 

retailer  productions  closest  but  prior  to  — , — , . . . , and  

no  no  no 

and  then  compute  the  average  echelon  inventory  as  if  these  were  the 
warehouse  production  points. 
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This  bound  can  be  achieved  in  certain  instances.  Consider  the  multiple 

cycle  policy  with  (n^n^n,,)  = (3,4,5).  The  warehouse  produces  at  either 
T T T 

— or  — . If  ■>  Q2>  the  warehouse  produces  at  — . Also,  the  ware- 

T 3 

house  produces  at  either  — or  T(— ).  If  Q2  >_  Q , the  warehouse  produces 

3 T 3 

at  T(— ).  Thus  if  Q„  = Q, , the  warehouse  produces  at  — and  T(— ) and 
5 2 1 4 5 

the  lower  bound  described  above  is  achieved. 

If  Q1  = Q2  in  the  multiple  cycle  policy  (3,4,5)  then  the  average 

1 3 

warehouse  echelon  inventory  is  Qq(—  + 20^ • However,  for  this  case 
L(3,4,5)  = Q0*min(i  + i,  + i)  = Qq(^-  + j).  Thus  L(nQ .n^ . . . ,nN)  may 
not  be  a lower  bound  for  general  multiple  cycle  policies  although  it  is 
an  easily  computed  approximation. 

The  above  results  give  us  a way  to  express  the  average  warehouse 
echelon  inventory,  or  a lower  bound  on  that  quantity,  for  multiple  cycle 
policies.  Given  such  an  expression,  we  can  write  out  a lower  bound  on  the 
average  cost  per  unit  time  for  a multiple  cycle  policy  solely  as  a function 
of  the  policy  (nQ .n^ . . . ,n  ) . We  can  then  compute  the  best  T*(nQ ,n  , . . . .n^ 
and  the  cost  C*(nQ ,n^ , . . . ,nN)  as  was  done  earlier  for  the  multiple  cycle 
policy  (2,2,3)  in  equations  (2.52)  and  (2.53).  We  can  then  compute  the 
best  unrestricted  n^'s  and  perform  the  same  iterative  procedure  as 
was  discussed  earlier  for  the  determination  of  single  cycle  policies.  The 
only  difference  is  that  we  must  ensure  n^  _>  nQ  for  all  i. 

Using  this  iterative  procedure , 50  of  the  randomly  generated  problems 
used  to  test  the  heuristic  procedures  were  evaluated  to  find  multiple  cycle 
policies  with  lower  average  cost  per  unit  time  than  the  best  single  cycle 
policy.  In  all  50  problems,  no  such  multiple  cycle  policy  was  found. 
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There  may  indeed  be  such  policies  but  none  were  found  with  this  heuristic 
search  technique. 

2 . 13  Summary 

There  appear  to  be  few  cases  where  one  can  guarantee  the  overall 
optimality  of  single  cycle  policies.  We  have  demonstrated  that  single 
cycle  policies  need  not  be  optimal  among  the  class  of: 

(1)  all  policies; 

(2)  separate  retailing  and  single  cycle  policies; 

(3)  stationary  policies; 

(4)  stationary  policies  satisfying  properties  P1-P5; 
and  (5)  stationary  renewal  policies. 

However,  single  cycle  policies  are  optimal  for  the  classes  of  problems 
where : 

(1)  all  retailers  are  identical; 

(2)  all  retailers  have  identical  ratios; 

and  (3)  the  overall  optimal  policy  is  stationary. 

In  many  of  the  proofs  presented,  non-stationary  policies  were  con- 
structed with  lower  cost  than  single  cycle  policies.  But,  it  is 
precisely  the  stationarity  of  single  cycle  policies  that  makes  them 
appealing  when  one  considers  the  practical  concerns  of  implementation. 

The  single  cycle  policy  is  not  only  stationary  in  lot  size  but  in  time  as 
well.  Each  facility  not  only  produces  the  same  lot  size  each  time  it 
produces  but  the  time  between  these  productions  is  constant  as  well. 

This  property  is  not  present  in  multiple  cycle  policies.  Also,  the 
presence  of  renewal  points  in  cycling  policies  offer  the  system  opportunities 
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to  restart  the  system  with  a revised  policy  should  any  of  the  system 
parameters  change.  Furthermore,  the  optimal  single  cycle  policy  or  near 
optimal  single  cycle  policy  is  easily  computed.  Graves  and  Schwarz  (1977) 
have  presented  a method  for  determining  good  single  cycle  policies  in 
general  arborescent  production  inventory  by  computing  single  cycle 
policies  for  many  one  warehouse,  N retailer  systems.  In  a procedure 
such  as  this,  the  efficient  computation  of  good  single  cycle  policies  is 
crucial. 

The  empirical  evidence  also  suggests  that  although  one  cannot  guarantee  i 

the  optimality  of  single  cycle  policies,  they  do  appear  to  be  close  to 

the  optimum,  at  least  for  small  problems.  Further  research  is  necessary  j 

to  establish  exact  bounds  on  the  difference  in  costs  between  the  best 

< 

single  cycle  solution  and  the  true  optimum  for  any  general  N retailer 
prob lem . 


¥ 


CHAPTER  3 

DYNAMIC  ONE  WAREHOUSE,  N RETAILER  PROBLEMS 

3.1  Introduction 

In  this  chapter  we  examine  the  dynamic  one  warehouse , N retailer 
production  inventory  system.  We  begin  by  introducing  the  basic  assump- 
tions and  notation.  We  review  the  Wagner-Whitin  and  Zangwill  algorithms 
for  the  single  facility  and  one  warehouse,  one  retailer  problems, 
respectively.  We  show  how  to  represent  the  one  warehouse,  N retailer 
problem  as  a single  source,  concave  cost,  network  flow  model.  A dynamic 
programming  algorithm  is  presented  for  this  case.  Easily  computed  upper 
and  lower  bounds  for  the  optimal  solution  value  are  developed  and  several 
special  cases  are  discussed. 

3 . 2 Basic  Model 

The  basic  structure  of  facilities  is  identical  to  that  discussed  in 
Chapter  2.  The  warehouse  will  be  referred  to  as  facility  0 and  the  N 
retailers  will  again  be  referred  to  as  facilities  1 through  N, 
respectively. 

In  the  dynamic  demand  case,  the  demand  at  retailer  i is  assumed  to 
occur  at  a known  rate  r^  in  period  t,  where  t = 1,2,...,T.  The  demand 
rate  may  differ  from  one  facility  to  another  and  from  one  time  period  to 
another.  We  will  assume  that  the  warehouse  experiences  no  external  demand. 

Inventory  is  maintained  between  stages  and  the  output  of  the  warehouse 
in  period  t can  be  used  by  any  of  the  retailers  in  the  same  period.  Time 
lags  can  be  included  by  redefining  the  decision  variables  to  account  for 

74 


75 


the  lead  time.  Also,  we  assume  that  backorders  are  not  allowed.  At  each 
facility,  chore  is  a production  cost  and  a holding  cost  which  may  differ 
from  one  period  to  the  next.  In  particular,  let 

x^  = Production  at  facility  j in  period  t, 

= On-hand  inventory  at  facility  j at  the  end  of  period  t, 

P^(x^)  = Production  cost  for  facility  j in  period  t, 

H^(I^)  = Holding  cost  for  facility  j in  period  t,  where 
j = 0,1,. . . ,N  and  t = 0,1,2,. ..,T. 

In  order  to  be  consistent  with  the  previous  research  on  related 
problems,  we  charge  the  holding  cost  at  the  warehouse  against  the  on-hand 
inventory  and  not  the  echelon  inventory  as  in  the  continuous  demand  case. 
We  will  also  assume  that  the  production  costs  and  holding  costs,  P^(x^) 
and  H^(I^)  are  concave  and  non-decreasing  for  non-negative  arguments. 
Without  loss  of  generality,  we  can  also  assume  that  1^  = 1^,  = 0 for 
all  j . 

The  objective  function  is  the  sum  of  the  holding  and  set-up  costs. 
The  constraints  reflect  the  inventory  balance  and  the  no  backlogging 
assumption.  We  can  therefore  express  the  problem  of  determining  the 
best  sequence  of  production  and  inventory  quantities  by: 
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T N . . . . 

Minimize  z = £ £ [P^(x^)  + H^(I^)]  (3.1) 

t=l  j=0  r T r * 

N 

0 0 0 r i 

subject  to  It  = It_1  + x - 2,  x^>  t = (3.2) 

j-1 

ij  = lj-i  + xj  " rt’  ^ (3.3) 

t = 1,...,T 

xj  > 0,  ij  > 0 (3.4) 

We  will  derive  a dynamic  programming  algorithm  for  this  problem  which 
is  a generalization  of  Zangwill's  (1963)  algorithm  for  the  serial  produc- 
tion inventory  system.  In  addition,  we  will  develop  uDper  and  lower  bounds 
on  the  optimal  solution  value  for  this  problem  which  made  use  of  Wagner 
and  Whitin' s (1958)  solution  technique  for  the  single  facility  problem. 
Before  discussing  the  algorithm  for  the  one  warehouse,  N retailer  system, 
it  is  therefore  necessary  to  review  these  two  previously  examined  cases. 

3 . 3 A Dynamic  Programming  Algorithm  for  the  Single  Facility  Problem 

In  this  section,  we  will  review  the  Wagner-Whitin  (1958)  solution 
technique  for  the  single  facility  problem.  The  algorithm  exploits  the 
structure  of  the  problem  when  determining  an  optimal  solution. 

For  the  single  facility  problem,  we  may  rewrite  the  mathematical 
programming  formulation,  (3.1)-(3.4),  as 
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T 

z = [ [P*(xb  + H^dh]  (3.5) 

t=l 

ij  = lj_1  + xj  - r£,  t = 1 , . . . ,T  (3.6) 

xj  >_  0,  ij  > 0,  t = 1,..  . ,T . (3.7) 

Since  the  objective  function  is  concave,  the  solution  to  the  above 
problem  must  occur  at  an  extreme  point  of  the  set  of  feasible  solutions. 

In  this  case, 

xj-lj  ]_  = 0 for>  all  t i-n  an  optimal  solution.  (3.8) 

Therefore,  x^  e {0,r^,r^+r^+1,. . . ,r^+r^+1+. . .+r^,}.  The  dynamic  programming 
algorithm  to  be  developed  will  take  advantage  of  this  form  for  an  optimal 
solution. 

Let  f (y)  = minimum  cost  for  the  first  k periods  when  the 

on-hand  inventory  at  the  end  of  period  k is  y. 

Then 


Minimize 


subject  to 


fk(y)  = min(Pk(xk)  + H^(y)  + > • (3-9) 


Since  x^*!^  ^ = 0, 


fk(y)  = min( 


+ + fk-i(o)> 


Ip^(O)  ♦ HjJCy)  + 


(3.10) 


Furthermore 
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fkty)  = min  Y^d)  where  (3.11) 

i 


Thus  Y^d)  = the  cos1:  Peri°ds  1 through  k under  the 

assumptions  that  an  optimal  policy  is  followed 
for  the  first  i-1  periods  and  the  inventory 
at  the  end  of  period  i-1  is  zero,  and  an 
order  arrives  at  the  beginning  of  period  i 
satisfying  the  demand  for  periods  i through  k 
with  the  on-hand  inventory  at  the  end  of  period 
k equalling  y . 

We  need  not  compute  f^  other  than  for  y = 0.  Thus  the  task  of  finding 
the  optimal  solution  requires  less  effort  than  applying  a standard  back- 
ward recursion  dynamic  programming  algorithm.  However,  for  the  one 
warehouse,  one  retailer  problem  the  backward  recursion  Zangwill  algorithm 
(.1969)  will  be  developed  in  the  next  section. 

3.4  A Dynamic  Programming  Algorithm  for  the  One  Warehouse,  One  Retailer 
Problem 

In  this  section,  we  review  a dynamic  programming  algorithm  for  the 
one  warehouse,  one  retailer  production  inventory  system.  Zangwill  (1969) 
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has  presented  a dynamic  programming  algorithm  for  the  serial  deterministic 
demand  production  inventory  system.  Clearly,  if  the  serial  system  has 
only  two  stages,  then  this  system  is  identical  to  the  one  warehouse, 
one  retailer  problem  described  in  the  preceding  section.  Zangwill  has 
shown  that  the  serial  system  may  be  represented  by  a single  source 
network  flow  model.  To  illustrate  this,  consider  Figure  3.1  which 
depicts  the  single  source  network  where  T = 4 and  node  (j,t) 
represents  facility  j and  period  t. 
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Clearly,  the  conservation  of  flows  in  the  network  in  Figure  3.1  ensures 
that  the  constraints  (3.2)  and  (3.3)  are  satisfied.  By  charging  the 
appropriate  cost  for  the  flow  along  each  of  the  arcs , it  is  clear  that 
finding  the  minimum  cost  solution  to  the  mathematical  programming  problem 
(3.1)-(3.4)  with  N = 1,  is  equivalent  to  finding  the  minimum  cost  flow 
in  the  network  representation. 

Since  the  cost  functions  are  concave  on  the  arc  flows  and  the  network 
has  a single  source,  we  know  that  at  least  one  "extreme  flow"  is  an 
optimal  solution  (Zangwill,  1968).  That  is,  at  most  one  arc  entering 
a node  will  have  a positive  flow. 

Therefore,  the  flow  into  each  node,  f,  can  be  expressed  by 

b 

f = 7 rT  where  t<a<b<T.  (3.13) 

.‘*1  — — — 
i=a 


In  this  case,  it  is  said  that  the  extreme  flow  has  the  "exact  requirements" 
property.  Clearly,  for  any  node  (l,t),  a = t in  (3.13)  to  prevent 
backorders.  Thus, 


(3.14) 


r 1 

From  (3.13),  if  the  flow  into  (0,t)  is  l r.,  tl.a^.b<_T, 

i=a  1 

then  there  are  two  cases  to  be  considered:  a = t and  a > t.  If 
a = t,  then  x^  > 0 to  insure  that  r^  is  satisfied.  Thus,  if  a = t, 
there  exists  an  integer  c,  t < c <_  b such  that 


4* 
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1 V 1 , t°  V 1 

Xt  = / ri’  and  Xt  = . i V 


(3.15) 


i=c+l 


If  a > t,  then  = 0 or  the  two  arcs  into  node  (l,t)  would  have 
positive  flow  and  therefore 


of1 

:t  = / V 


(3.16) 


Define  K^(a,b)  to  be  the  optimal  cost  of  completing  the  production 

b ! 

schedule  if  l r.  units  enter  node  (j,t).  Also  let 
i=a  1 


Pj(a,b)  = ?3t  l 
Hj(a,b)  = I r^  . 


If  a = t in  (3.13),  for  the  retailer  node  (l,t),  then  by  (3.14), 


K*(t,b)  = Hd(t+l,b)  + K^+1(t+l,b) 


l H^(i+l,b) , 
i=t  1 


(3.17) 


which  holds  for  1 <_  t <_  T-l,  t^_b^_T;  if  t = T then  K^,(T,T)  = 0. 
If  a = t in  (3.13),  for  the  warehouse  node  (0,t),  then  by  (3.15), 


K°(t,b)  = min  {pj(t,c)  + H°(c+l,b)  + K*(t,c)  + K°  (c+l,b)>,  (3.18) 

t t<c<b  T r T r 1 

for  1 < t < T-l  and  t<_b<_T;  if  t = T,  kJ(T,T)  = pJ(T,T). 
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If  a > t in  (3.13)  for  the  warehouse  node  (0,t),  then  by  (3.16), 

K°(a,b)  = H°(a,b)  + K°+1(a,b),  (3.19) 

for  2 <_  t <_  T-l  and  t < a <_  b <_  T.  If  t = l,  then  a = 1 and 

K^(l,b)  from  (3.18)  can  be  applied. 

To  complete  the  algorithm,  we  must  determine  the  decisions  at  the 

source  node,  i.e.  how  much  to  produce  at  the  warehouse  in  each  period. 

T 

S — i 

Let  Kt(a,T)  = minimum  cost  of  shipping  £ r units  from  the  source 

t=a 

to  satisfy  the  requirements  at  nodes  (l,a) ,(l,a+l) ,. . . ,(1,T) , that 
is,  the  optimal  cost  of  satisfying  demands  r^,...,rj  through  production 
in  periods  t,t+l,...,T.  Demand  is  assumed  satisfied  for  the  first  a-1 
periods  from  production  in  the  periods  prior  to  t.  Also  if  a = t,  we 
must  produce  in  period  t at  the  warehouse. 

Therefore , 

K^(a,T)  = minimum  (P°(a,c)  + K°(a,c)  + (c+l,T)},  (3.20) 

max(t ,a-l)£c£T  + 

where  2 <_  t <_  T-l  and  t <_  a <_  T.  if  c = a-1,  x°  = 0 by  convention 

and  if  t = 1,  we  use  (3.20)  only  with  a = 1.  . 

We  begin  the  algorithm  with  K^(T,T)  = 0 and  work  backward  until 
s s 

K^(1,T)  is  evaluated.  K^(1,T)  is  the  minimum  cost  attainable  in  this 
network. 

Zangwill  (1966,1969)  has  extended  the  single  facility  model 
discussed  in  Section  3.3  and  the  one  warehouse,  one  retailer  model 
discussed  above  to  include  the  possibility  of  backlogging  of  external 
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i 


demand.  Each  retailer  node  (1,T)  in  Figure  3.1  is  revised  to  include 
additional  arcs  to  represent  the  flow  back  in  time  (Figure  3.2).  If  we 
let 


.1+ 


I = inventory  on-hand  at  the  end  of  period  t at  retailer  1, 


,1- 


inventory  short  at  the  end  of  period  t at  retailer  1, 


H^+(I^+)  = inventory  carrying  cost  associated  with  I^+, 
and  (1^  ) = shortage  cost  associated  with 


© 


I1  + 
HH 


'It 


il+ 


FIGURE  3.2  NODE  (1,  t)  WITH  BACKLOGGING 

We  will  omit  the  development  of  the  recursion  because  we  will  only 
extend  the  no  backlogging  case  to  the  one  warehouse,  N retailer  system 
An  extension  to  include  backlogging  is  possible.  But  as  we  shall  see  in 
the  next  section,  the  one  warehouse,  N retailer  algorithm  is  quite 
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complicated  computationally  and  such  extensions  will  only  further  exacer- 
bate this  problem. 

3 . 5 A Dynamic  Programming  Algorithm  for  the  One  Warehouse , N Retailer 
Problem 

It  is  possible  to  generalize  the  one  warehouse,  one  retailer  dynamic 
programming  algorithm  presented  in  Section  3.4  to  solve  the  one  warehouse, 
N retailer  system.  Veinott  (.1969)  first  noted  that  this  generalization 
was  possible  but  he  failed  to  point  out  the  basic  differences  between 
the  generalized  algorithm  and  the  Zangwill  procedure.  Veinott  also 
presented  an  efficient  algorithm  for  the  general  arborescent  problem 
under  a severe  set  of  restrictions  on  the  cost  structure  that  is  similar 
to  Love's  nested  algorithm  (1972)  for  facilities  in  series. 

We  begin  by  representing  the  one  warehouse,  N retailer  system 
by  a single  source  network  flow  model.  To  illustrate  this,  consider 
Figure  3.3  which  depicts  the  single  source  network  for  two  retailers 
where  T = 4 and  node  (j,t)  represents  facility  j and  period  t. 

Clearly,  the  conservation  of  flows  in  the  network  in  Figure  3.3 
ensures  that  the  constraints  (3.2)  and  (3.3)  are  satisfied.  By  charging 
the  appropriate  cost  for  the  flow  along  each  of  the  arcs , it  is  clear 
that  finding  the  minimum  cost  solution  to  the  mathematical  programming 
problem  (3.1)- (3. 4)  is  equivalent  to  finding  the  minimum  cost  flow  in 
the  network  representation. 

Since  the  cost  functions  are  concave  on  the  arc  flows  and  the  network 
has  a single  source , we  know  that  at  least  one  extreme  flow  is  an  optimal 
solution.  That  is,  at  most  one  arc  entering  a node  will  have  a positive 


flow. 
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b 

f = 7 where  t < b < T , 

i=t  1 


and  therefore 


I]  = 


i=t+l 


r^ 

1 


t < b < T. 


b = t 


(3.21) 


Define  r(t,b),  as  before,  to  be  the  optimal  cost  of  shipping 
b . Z 

£ r?  units  from  node  (j,t)  to  their  respective  destinations.  Then, 
i=t  1 

it  is  clear  that 


Kj(t,b)  = Hj(t+l,b)  + Kj+1(t+l,b)  , 


(3.22) 


.which  holds  for  j > 0,  1 <_  t <_  T-l,  t^_b<_T;  if  t = T,  then 

K^,(T ,T ) = 0 for  j > 0.  Also, 

Kj(t,t)  = Hj(t+l,t)  + Kj+1(t+l,t),  (3.23) 

which  by  convention. 


= 0 + 0. 


The  recursion  is  somewhat  more  complicated  for  the  warehouse  nodes 
(0,t).  The  extreme  flow  property  demonstrates  that  the  flow,  f,  into 
a node  (0,t)  can  be  expressed  by 


b 

*■  l 


ri  + 

i=a  i=c 


r 2 

.1  ri’ 


(3.24) 
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where  there  are  3 cases  for  the  values  of  a,  b,  c and  d.  Either 
(i)  * .1  a b .1  T and  t _<  c < d < T, 

or  (ii)  a = b+1,  b <_  T and  t <_  c <_  d < T, 

or  (iii)  t _<  a <_  b _<  T and  c = d+1,  d _<  T. 

Case  (i)  is  exactly  analogous  to  the  flow  into  a warehouse  node  for  the 

b 

one  retailer  case.  By  convention,  £ r^  = 0,  for  j = 1,2,  and 

i=b+l  1 

case  (ii)  and  (iii)  are  necessary  because  the  flow  into  a warehouse  node 

may  only  consist  of  the  sum  of  requirements  for  a single  retailer. 

0 b 1 d 2 
Let  Kt(a,b,c,d)  be  the  optimal  cost  of  shipping  £ r.  + £ rf 

i=a  1 i=c  1 

units  to  their  respective  destinations.  Also  define 


H (a,b,c,d)  = H ( 7 rd  + 7 r.)  and 

t t . L l . u l 

i=a  i=c 

P°(a,b,c,d)  = pJ(  l rd  + [ r?). 

i=a  i=c 


By  convention, 


H°(b+l,b,d+l,d)  = P°(b+l,b,d+l,d)  = 0. 


We  begin  by  considering  case  (i)  t ^ a b <_  T,  t < d < T 
for  which  there  are  four  subcases,  since  both  a and  c either  equal 
t or  are  greater  than  t. 


(i-1)  If  a = t and  c = t,  there  exist  integers  f and  g,  t <_  f <_  b, 
t < g < d such  that 


„1  r 1 v2  ? 2 . _0  r 1 r 2 

xt  = L ri,  xt  = I ri  and  I = I r.  + I r 

i=a  i=c  i=f +1  i=g+l 


(3.25) 
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( i—  2 ) If  a = t and  c > t,  there  exists  an  integer  f , t < f < b, 
such  that 


1 V 

xt  3 J 


i=a 


1 

ri’ 


0 and 


t°  r 1 r 2 

:t  = L ri  + I ri- 

i=f+l  i=c  1 


(3.26) 


(i-3)  If  a > t and  c = t,  there  exists  an  integer  g,  t < g < d, 

such  that 

x^  = 0,  x^  = 7 r?  and  1^  = 7 r^  + 7 r?.  (3.27) 

t t.Li  t .**  l . L , l 

i=c  i=a  i=g+l 

(i-4)  If  a > t and  c > t,  then 

i 9 n b i d o 

x^  = 0,  x = 0 and  I = 7 rT  + 7 rT.  (3.28) 

t t t .L  l .L  l 

i=a  i=c 


We  can  now  determine  the  recursive  equations  for  each  of  these  4 
subcases  of  case  (i). 

(i-1)  If  a = t and  c = t,  by  (3.25), 


K°Ct,b,t,d)  = min  fp*(t,f)  + K^(t,f)  + H°(f+l,b,g+l,d) 

t<f<b  ( T >;(3.29) 


t-g-d  L + P^(t,g)  + K^(t,g)  + K°+1(f+l,b,g+l,d)_ 


for  1 <_  t <_  T-l; 

if  t = T,  k£(T,T,T,T)  = P*(T,T)  + P*(T,T>. 
( i—  2 ) If  a = t and  c > t,  by  (3.26), 
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K°(t,b,c,d)  = min  jpht.f)  + K^(t,f ) + H°(f+l,b,c,d) 


t<f  <b 


for  1 < t < T-l. 


+ Kt+1(f+l,b,c,d) 


(i-3)  If  a > t and  c = t,  by  (3.27), 


, (3.30) 


K°(a,b,t,d)  = min  fp^(t,g)  + K^(t,g)  + H®(a,b,g+l,d)! 

t<g<d  ( r x T /,  (3.31) 


Kt+1(a,b,g+l,d) 


for  1 < t < T-l. 


(i-4)  If  a > t and  c > t,  by  (3.28), 


K°(a,b,c,d)  = H°(a,b,c,d)  + K°+1(a,b,c,d). 


(3.32) 


Thus  far  the  N retailer  case  is  a straightforward  extension  of 
Zangwill's  one  retailer  algorithm.  But,  the  development  of  recursion 
equations  for  cases  (ii)  and  (iii)  is  somewhat  more  complicated.  Let  us 
examine  case  (ii)  more  carefully. 

b x 

Case  (ii)  corresponds  to  £ r.  = 0 units  for  retailer  1 and 
d 2 i=b+l  1 

£ r.  units  for  retailer  2 carried  into  the  warehouse  node  (0,t). 
i=c 

There  is  only  one  possible  way  to  arrive  at  such  a situation.  In  either 
the  recursion  (3.29)  or  (3.30),  for  some  preceding  time  period,  f was 
chosen  equal  to  b.  That  is,  the  requirements  up  to  and  including  period 
b were  produced  and  then  held  at  retailer  1.  Clearly,  the  retailer  1 
requirements  for  period  b+1  and  after  have  not  yet  been  produced  at 


the  warehouse. 
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2 

Let  us  suppose  now  that  b = t-1.  Then  in  case  (ii),  £ r.  units 

i=c 

flow  into  node  (0,t)  but  as  discussed  above  the  requirements  for 
retailer  1 for  period  t and  after  have  not  yet  been  produced. 

Clearly,  to  maintain  feasibility  at  least  r^  units  must  be  produced 
at  the  warehouse  in  period  t.  This  implies  that  the  two  arcs  into  node 
(0,t)  both  have  positive  flow  which  is  impossible  by  the  property  of 
extreme  flows.  A similar  argument  for  case  (iii)  illustrates  that 
d = t-1  yields  the  same  contradiction.  To  ensure  that  the  dynamic 
programming  algorithm  only  examines  extreme  flows  we  define 


K^(t,t-l,c,d)  = <*>  for  c ^ d and 


K (a,b,t,t-l)  = for  a ^ b. 


(3.33) 


Given  these  definitions,  we  may  now  write  the  recursive  equations 

for  cases  (ii)  and  (iii).  For  case  (ii)  there  are  two  subcases  since 

2 

either  c = t or  c>t.  If  c = t,  x^  > 0 and  we  may  write 


K°(b+l,b,t,d)  = min  jV(t,g)  + K^(t,g)  + K°  (b+l,b,g+l,d)) 
t<g<d  < x r 

l + H°(b+l,b,t,g) 


(3.34) 


If  c > t,  x^  = 0 and  we  may  write 


K°(b+l,b,c,d)  = H°(b+l,b,c,d)  + K°+1(b+l,b,c,d) . 


(3.35) 


Similarly,  for  case  (iii)  the  two  subcases  are  whether  a = t or  a > t. 


If  a = t,  x^  > 0 and  we  may  write 
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K°(t,b,d+l,d)  = min  fp*(t,f)  + K*(t,f)  + K°  (f+l,b,d+l,d) 
tj<f<b/  T 

\ + H°Ct,f,d+l,d) 


. (3.36) 


If  a > t,  = 0 and  we  may  write 


K°(a,b,d+l,d)  = H°(a,b,d+l,d)  + K^+1(a,b,d+l,d). 


(3.37) 


Finally, 


K^(b+l,b  ,d+l,d)  = 0. 


(3.38) 


To  complete  the  algorithm,  we  must  determine  the  decisions  at  the 

source  node,  s,  i.e.  how  much  to  produce  at  the  warehouse  in  each 

T T 

s * 1 * 2 

period.  Let  K (a,T,c,T)  = minimum  cost  of  shipping  £ r + £ r^ 

t=a  t=c  C 

units  from  the  source  to  satisfy  the  requirements  at 

(l,a) ,(l,a+l) , . . . ,(1,T) ,(2 ,c) ,(2 ,c+l) , . . . ,(2 ,T) , that  is,  the  optimal 

1 12  2 

cost  of  satisfying  demands  r ,...,r  ,r  ,...,r  through  productions 

3 I C 1 

in  periods  t,t+l,...,T.  Demand  for  each  retailer  is  assumed  satisfied 
from  production  in  the  periods  prior  to  t.  Also  if  a = t or  c = t, 
we  must  produce  in  t at  the  warehouse. 

Therefore 


K®(a,T,c,T)  = minimum  |P°(a,f,c,g)  + K^(a,f,c,g) 
max(t,a-l)<f<T  / 


raax(t,c-l)<g<T  ^ + ^(f+l.T.gtl.T) 


(3.39) 


where  2 <_  t <_  T-l  and  t<^a,c^_T.  If  t = 1,  we  use  (3.39)  only  with 
a = c = 1. 


i 
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1 2 

We  begin  the  algorithm  with  K^,(T,T)  = KT(T,T)  = 0 and  work  backwards 
for  each  of  the  retailers , followed  by  the  warehouse  and  then  the  source 

g 

node  until  we  compute  K^(1,T,1,T)  which  is  the  minimum  cost  attainable 
in  the  network. 

Clearly,  the  computational  effort  involved  in  computing  the  minimum 
cost  attainable  in  the  network  increases  exponentially  in  the  number  of 
retailers.  However,  the  properties  of  an  extreme  flow  allow  us  to  cut 
down  the  computational  effort  considerably. 

Theorem  3.1.  If  the  flow  into  a node  (0,t)  is  an  extreme  flow  and  is 
b i d 2 

equal  to  £ r . + £ r . , where  t <_  a <_  b <_  T and  t <_  c <_  d <_  T , 

i=a  i=c 

then  (a,b)  intersects  (c,d)  on  the  time  axis,  i.e.  it  is  impossible 
for  either  a<_b<c<_d  or  c<_d<a<_b  to  hold. 

Proof  of  Theorem  3.1. 

d 2 

Assume  a <_  b < c <_  d.  Clearly,  the  \ rf  units  must  be  carried 

i=c  1 

at  the  warehouse  level  until  period  c.  This  must  be  the  case  since  to 

2 

insure  feasibility,  I >0  for  t-1  t <_  c-1  and  for  any  of  the 

d 2 T 

\ rf  units  to  be  carried  at  the  retailer  level  prior  to  period  c would 

i=c  2 

require  x^  > 0 for  some  t <_  t <_  c-1.  This  is  impossible  by  the  property 

of  extreme  flows.  Therefore,  1°  > 0 for  t,  where  t < t < c-1. 

T — — 

By  assumption  b < c and  must  be  satisfied  by  period  b+1 

to  maintain  feasibility.  Therefore,  x^  > 0 for  t,  where 
t < t < b+1  < c.  Thus  x0.*I°,  , > 0 for  some  t',  where  t+1  < t*  < c 
which  contradicts  the  extreme  flow  property.  Therefore,  it  cannot  be 
the  case  that  a <_  b < c <_  d.  The  proof  for  the  case  where  c d < a <_  b 
simply  reverses  the  role  of  the  retailers  in  the  above  argument.  □ 
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It  may  be  possible  to  improve  the  computational  efficiency  of  the 
dynamic  programming  algorithm  further  through  the  use  of  bounds  on  the 
optimal  solution  value.  Such  bounds  will  be  discussed  in  the  following 
section. 

3.6  Bounds  on  the  Optimal  Solution  Value  for  One  Warehouse,  N Retailer 
Problems 

The  dynamic  programming  algorithm  presented  in  Section  3 . 5 computes , 
for  each  node  (j,t),  the  optimal  cost  of  shipping  a group  of  exact 
requirements  from  that  node  to  their  ultimate  destinations.  If  the 
cost  of  such  a partially  completed  production  inventory  plan  exceeds  the 
cost  of  a known  feasible  solution  to  the  entire  problem,  then  we  can 
omit  consideration  of  any  completions  of  this  partial  plan.  This  is 
clearly  the  case  if  all  the  cost  functions  are  non-negative,  which  is  a 
realistic  assumption. 

Obviously,  the  closer  the  known  feasible  solution  value  is  to  the 
optimum,  the  more  likely  it  is  that  this  bounding  procedure  will  enable 
us  to  eliminate  prospective  solutions  and  thus  save  computational  effort . 
We  will  discuss  two  methods  of  determining  feasible  solutions  that  will 
hopefully  have  this  property.  To  simplify  the  discussion,  we  will  again 
only  explicitly  consider  the  class  where  N = 2.  The  discussion  is 
easily  extended  for  larger  values  of  N. 

The  first  method  consists  of  N+l  = 3 applications  of  the  Wagner- 
Whitin  single  facility  algorithm  discussed  in  Section  3.3.  Begin  by 
solving  for  the  best  production  inventory  plan  for  the  retailers,  where 
each  one  is  treated  as  an  independent  single  facility  problem.  Let  us 


now  assume  that  this  portion  of  the  one  warehouse , 2 retailer  plan  is 

fixed.  Consider  the  warehouse  node  (0,t).  Define  = £ x^  which 

j=l 

is  determined  by  the  choice  of  the  production  plans  for  the  retailers. 
Therefore,  given  the  production  plans  for  the  retailers,  the  problem 
reduces  to  a single  facility  problem  with  requirements  r^  and  costs 
and  P®  which  can  again  be  solved  by  the  Wagner-Whitin  algorithm. 
This  plan  guarantees  feasibility  by  satisfying  all  the  demands  at 
the  retailers  and  insuring  that  the  inventory  balance  constraints  at  the 
warehouse  nodes  are  met.  It  is  also  clear  that  this  feasible  solution 
will  be  an  extreme  flow.  Furthermore,  the  value  of  any  feasible  solution 
is  an  upper  bound  on  the  optimal  solution  value.  This  method  of  decom- 
posing the  system  is  similar  to  the  method  proposed  by  Kalymon  (1972) 
for  general  arborescent  systems. 

Since  the  optimization  at  the  retailer  level  does  not  take  the 
warehouse  costs  into  account , the  feasible  solution  obtained  by  this 
technique  need  not  be  the  overall  optimal  solution.  We  can,  for  a 
particular  cost  structure  at  the  warehouse,  revise  the  production  costs 
at  the  retailer  to  reflect,  at  least  in  part,  the  warehouse  costs. 
Consider  the  node  (0,t)  in  the  single  source,  network  flow  representa- 
tion of  the  one  warehouse,  N retailer  problem  (Figure  3.4). 

If  Cx^+x^)  > 0,  then 


'At* 


+ X 


2 

t 


to  maintain  feasibility. 


By  the  property  of  extreme  flows 


0 


and  so  either 


Therefore,  a cost  of  at  least  min(Ht  1(xt+xt^»  Pt(*t+xt))  is  incurred 

1 2 

at  the  warehouse  level  by  fixing  + x at  the  retailer  level. 


FIGURE  3.4  NODE  (0,0  IN  THE  NETWORK  FLOW 
REPRESENTATION 


if 


H°Cx)  = 


0 if  x = 0, 

and 

H°  + h°x  if  x > 0, 


r 0 if  x ■ 

l Kj  ♦ pjx  if  x > 0, 


where  t = 1 


T,  then  we  can  revise  the  production  costs  at  the 
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retailers  to  reflect  the  cost  that  must  be  incurred  at  the  warehouse  level 
by  fixing  the  production  levels  at  the  retailers.  That  is,  set 


if  x = 0 


K^(x)  = 


H°  K° 


for  i = 1,2.  (3.40) 


K^(x)  + min(-^-,-^-)  + min(h°x,p°x) , if  x 


The  H^/2  and  K^/2  are  present  to  insure  that  if  both  retailers  produce 
in  period  t,  only  a single  fixed  charge  of  either  or  is 

incurred  at  the  warehouse.  (A  similar  revision  in  costs  can  be  made  if 
the  warehouse  costs  have  either  a fixed  charge  at  zero  or  a linear  charge 
but  not  both.) 

Thus  we  may  summarize  the  technique  to  find  a feasible  solution  as 
follows : 

(1)  For  each  retailer  i,  use  the  Wagner-Whitin  algorithm 

with  costs  and  and  requirements  r^  to  compute 

a production  plan  for  each  retailer. 

0 12 

(2)  Compute  r^  = x + x^_  for  all  t. 

(3)  Solve  for  the  best  production  plan  at  the  warehouse 
using  a Wagner-Whitin  algorithm  with  costs  K°,  H®  and 
requ irement  s r ® . 

(4)  Compute  the  overall  cost  of  this  production  plan  using 
the  original  costs  for  the  network. 

An  alternative  technique  to  find  a feasible  solution  for  this 
system  is  analogous  to  the  separate  retailing  policy  for  the  continuous 
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demand  case.  Simply  split  the  one  warehouse,  N retailer  system  into  N 
one-retailer  systems  and  solve  each  of  these  systems  independently  using 
Zangwill's  algorithm.  This  technique,  however,  requires  much  more  compu- 
tational effort  than  the  method  described  above  which  used  the  Wagner- 
Whitin  procedure. 

We  can  also  determine  the  quality  of  the  feasible  solution  determined 
by  either  of  the  two  heuristics  described  above  by  comparing  the  value 
of  the  solution  to  the  value  of  a lower  bound  on  the  optimal  solution.  One 
possible  lower  bound  can  be  computed  by  summing  the  costs  of  the  Wagner- 
Whitin  plans  for  each  of  the  retailers,  omitting  the  warehouse  costs 
except  in  the  cost  K^.  That  is,  these  plans  are  determined  and  evaluated 
with  the  revised  production  cost  for  retailer  i.  Clearly,  no  plan 

can  be  found  with  revised  costs  lower  than  this  for  the  retailers  and 
was  chosen  to  reflect  the  minimum  cost  added  at  the  warehouse  for 
any  plan  at  the  retailers. 

An  alternative  lower  bound  may  be  determined  by  computing  the  Wagner- 
Whitin  solution  for  the  single  facility  with  production  costs  K°,  holding 
costs  H^Cx)  = minCH^Cx) ,H^(x) ,H^(x) ) and  requirements  r£  + r^.  We 
can  demonstrate  that  this  is  indeed  a lower  bound,  if  we  examine  the 
mathematical  programming  model  for  the  one  warehouse,  2 retailer  system, 
i.e.  (3.1)-(3.4)  with  N = 2.  If  we  define  = I?  + 1^  + I?,  we  can 
rewrite  the  inventory  balance  equations,  (3.2)  and  (3.3)  as 


x 


0 

t 


+ J 


0 

t-1 


_1  ♦ v,2 

rt  + V 


(3.41) 


By  the  definition  of  H,  it  is  clear  that 


Furthermore 


Ht(I°+lJ)  + H (I*)  by  concavity 


Similarly, 


H ( J . ) < H°(I°)  + Hj(lj)  + H*(I*). 
t 1 — t t tt  tt 


Therefore,  the  objective  function  value  for  this  relaxed  mathematical 
program  with  constraint  (3.41)  and  holding  cost  function  is  a 

lower  bound  on  the  optimal  solution  value  for  the  original  problem. 

These  lower  bounds  can  be  achieved  when  certain  restrictions  are 
placed  on  the  cost  structure  of  the  system.  These  special  cases  will 
be  discussed  in  the  following  section. 


3.7  Exact  Solutions  for  Restricted  Cost  Structures  for  One  Warehouse, 

N Retailer  Problems 

It  is  also  possible  to  use  the  Wagner-Whitin  algorithm  to  get  the 
optimal  solutions  for  certain  special  cases  of  the  cost  structure  in  the 
one  warehouse,  N retailer  system. 

Theorem  3.2.  If  K^(x)  = px,  then  the  optimal  solution  for  the  one 
warehouse,  N retailer  system  may  be  found  by  solving  N independent 
single  facility  problems,  one  for  each  retailer. 
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Proof  of  Theorem  3.2. 

If  K^(x)  = px,  we  may  ignore  the  per  unit  production  costs  altogether 

since  the  total  number  of  units  produced  at  the  warehouse  must  equal 
T N 

1 Jr.  Therefore,  the  cost  of  production  at  the  warehouse  may  be 
t=l  i=l  X 

assumed  to  be  0.  Clearly,  no  units  will  be  carried  at  the  warehouse 
from  period  t to  t+1  provided  is  non-negative.  Summing  the 

costs  for  the  N independent  Wagner-Whitin  solutions  is  clearly  a lower 

o N 
0 r 1 

bound  on  the  optimal  cost.  If  we  set  x = J x for  t = 1,2,...,T, 

X i=l  t 

then  the  additional  charge  for  this  part  of  the  production  plan  is  0 
and  the  overall  cost  of  the  policy  achieves  the  lower  bound.  □ 

Although  the  assumption  that  K^(x)  = px  may  seem  unrealistic, 
there  are  cases  where  such  a cost  structure  is  likely.  If  the  warehouse 
receives  deliveries  each  time  period  from  its  supplier,  then  there  may 
be  no  additional  fixed  cost  associated  with  placing  orders  in  a particular 
policy.  Similarly,  if  shipments  are  made  each  period  from  the  warehouse 
to  each  of  the  retailers,  then  there  may  be  no  additional  fixed  set-up 
cost  associated  with  a retailer  placing  an  order  in  a particular  policy. 

In  this  case,  we  may  show 

Theorem  3.3.  If  for  all  t and  i, 

K*(x)  = p*x,  and 
H°(x)  < H*(x), 

then  the  optimal  solution  for  the  one  warehouse,  N retailer  system  may 


Li 
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be  found  by  solving  the  single  facility  problem  with  costs  K°(x),  H°(x) 

0 N i t t 

and  requirements  r = £ r . 

X i=l 

Proof  of  Theorem  3.3. 

As  in  the  proof  of  Theorem  3.2  for  the  warehouse,  we  may  ignore 
the  production  costs  at  the  retailers  entirely  in  this  case.  Now 
suppose  that  I units  are  held  from  period  t to  t+1  at  retailer  i. 
Clearly,  we  can  decrease  the  preceding  retailer  production  by  I units. 
If  those  I units  are  carried  at  the  warehouse  until  period  t+1,  a 
lower  cost  is  incurred  since  H^(x)  <_  H^(x)  and  is  concave.  Adding 

a retailer  i production  in  period  t adds  nothing  to  the  cost. 
Therefore,  all  units  are  carried  at  the  warehouse  in  an  optimal  policy. 

This  case  is  analogous  to  the  second  lower  bound  on  the  optimal 
solution  value  discussed  in  the  preceding  section.  By  that  same  argu- 
ment, the  optimal  cost  for  the  single  facility  problem  with  costs 

. N 

K^(x),  H^(x)  = min(H^(x))  = H°(x),  and  requirements  r°  = £ r^  is 

i i=l 

a lower  bound  on  the  overall  optimal  cost.  Since  the  production  costs 

at  the  warehouse  are  ostensibly  zero,  this  lower  bound  is  achieved  by 
following  the  plan  for  the  warehouse  discussed  above  and  simply  setting 
x^  = r*  for  all  i and  t.  □ 

A similar  result  can  be  proven  if  retailer  i has  lower  holding 
cost  over  all  t than  the  warehouse  has. 

Corollary  3.4.  If  K°(x)  = K + p°x,  K^(x)  = p^x  for  all  t and  i 
and  for  each  i either 

(a)  H°(x)  <_  H*(x),  for  all  t and  x >_  0 

H^(x)  = h^*x  < h°*x  = H^(x),  for  all  t,  and  x ^ 0, 


or  (b) 
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then  the  optimal  solution  of  the  one  warehouse,  N retailer  system  may  be 
found  by  solving  an  appropriate  single  facility  problem. 

Proof  of  Corollary  3.4. 

As  in  Theorem  3.3,  we  may  ignore  the  retailer  production  costs  totally. 
If  H^(x)  <_  H^(x),  then  no  units  are  carried  at  retailer  i in  an  optimal 
solution.  If  for  all  i,  this  is  the  case  then  the  corollary  reduces 
to  Theorem  3.3. 

Suppose  H^Cx)  = h^x  < h°-x  = H°(x)  for  some  i and  suppose  r^ 
is  carried  at  the  warehouse  from  some  period  t to  some  period  x + 6. 
Clearly,  we  may  reduce  1^  by  r^  units  and  add  a retailer  i production 
in  period  x.  This  must  reduce  the  holding  costs  over  the  periods  x to 
x + 6 since  h^x  < h0<x.  Therefore,  no  units  are  carried  at  the  warehouse 
for  retailer  i. 

Thus,  the  requirements  for  a particular  retailer  are  either  carried 
exclusively  at  the  warehouse  or  exclusively  at  the  retailer  depending  on 
whether  case  (a)  or  (b)  holds  concerning  the  holding  costs.  Denote  the 
retailers  that  satisfy  (a)  by  1,2,..., a. 

It  is  also  clear  that  if  the  warehouse  produces  in  period  t,  then 
it  produces  for  all  the  retailers.  If  not,  then  the  holding  costs  may  be 
reduced  by  decreasing  the  preceding  warehouse  production  and  increasing 
the  production  in  period  t.  By  assumption,  this  cannot  increase  the 
production  costs  at  the  warehouse. 

If  each  time  the  warehouse  produces  it  produces  for  all  the  retailers, 
then  it  is  clear  that  the  system  has  common  regeneration  points.  There- 
fore, we  may  solve  a single  facility  problem  with  production  costs  of 

N i 

K°  and  requirements  of  £ r in  period  t,  provided  that  we  charge 
* i=l 
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the  holding  costs  properly.  In  particular,  if  the  requirements  for  period 

t+1  are  carried  from  period  t to  t+1,  the  system  holding  costs  are 
a . N . . 

H ( £ r^)  + £ hXT^.  By  suitable  redefinition  of  (3.12)  in  the 

i=l  i=a+l 

Wagner-Whitin  recursion,  we  can  solve  this  problem  by  that  algorithm  with 
some  additional  computation  but  no  increase  in  the  number  of  decisions.  □ 

As  was  mentioned  in  Chapter  1,  an  extension  of  the  Wagner-Whitin  model 
that  has  received  a great  deal  of  attention  is  the  problem  with  capacity 
constraints.  Florian  and  Klein  (1971),  Lambrecht  and  Vander  Eecken 
(1978)  and  Baker,  Dixon,  Magazine  and  Silver  (1978)  have  all  proposed 
algorithms  to  solve  this  particular  problem. 

Theorems  3.2,  3.3  and  Corollary  3.4  demonstrate  that  under  certain 
restrictions  on  the  cost  structure  that  one  can  solve  the  one  warehouse, 

N retailer  problem  by  the  solution  of  single  facility  problems  using  the 
Wagner-Whitin  algorithm.  Similarly,  if  there  are  production  capacities 
on  some  of  the  arcs,  then  for  certain  cost  structures  it  is  possible  to 
solve  the  one  warehouse,  N retailer  capacitated  system  by  the  solution  of 
single  facility,  capacity  constrained  problems. 

Corollary  3.5.  If  K°(x)  = px  and  there  are  capacity  constraints 
x^  c^»  i = 1* • • • »N,  t = 1, . . . ,T , then  the  optimal  solution  for  the 
one  warehouse,  N retailer  system  may  be  found  by  solving  the  N 
independent  single  facility,  capacity  constrained  problems. 

Proof  of  Corollary  3.5. 

As  in  the  proof  of  Theorem  3.2,  we  may  ignore  the  production  costs 
at  the  warehouse  entirely,  and  also  no  units  need  by  carried  at  the  ware- 
house. The  solution  for  the  general  one  warehouse,  N retailer  capacity 
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constrained  problem  is  bounded  below  by  the  sum  of  the  costs  of  the  optimal 
solutions  for  the  N independent  single  facility  capacity  constrained 
problem.  For  this  cost  structure  the  bound  is  achieved.  □ 

In  a similar  fashion,  we  may  demonstrate 

Corollary  3.6.  If  for  all  t and  i > 0 K^(x)  = p^x,  H^(x)  <_H^(x), 

and  x^  <_  c°,  then  the  optimal  solution  for  the  one  warehouse,  N 
retailer  capacity  constrained  system  may  be  found  by  solving  the 
capacity  constrained  single  facility  problem  with  costs  K^(x),  H^(x) 

0 N i 1 

and  requirements  r = £ r\  □ 

i=l 

A similar  restatement  of  Corollary  3.4  cannot  be  made  because  of 
the  need  to  shift  units  produced  at  the  warehouse  to  demonstrate  common 
regeneration  points.  The  capacity  constraints  on  the  production  arcs 
at  the  warehouse  preclude  this  argument. 

3.8  Summary 

In  this  chapter,  we  have  examined  the  dynamic  one  warehouse,  N 
retailer  production  inventory  system.  This  problem  is  a generalization 
of  the  2-stage  serial  system  discussed  by  Zangwill.  The  dynamic  program- 
ming algorithm  presented  extends  the  concepts  of  Zangwill' s algorithm 
to  this  more  general  problem.  There  are  some  computational  simplifica- 
tions that  can  be  made  and  upper  and  lower  bounds  for  the  optimal 
solution  value  are  discussed.  The  method  for  determining  upper  bounds 
is  a heuristic  technique  to  determine  "good"  feasible  solutions. 

Although  the  dynamic  programming  algorithm  may  seem  intractable 
for  problems  with  a large  number  of  retailers,  there  are  a number  of 
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special  cases  where  single  facility  algorithms  may  be  applied  to  determine 
the  optimal  solution.  Further  investigation  is  necessary  to  establish 
when  these  much  simpler  computational  techniques  yield  solutions  that 
are  either  optimal  or  sufficiently  close  to  optimal. 


CHAPTER  4 


CONCLUSION 

In  this  thesis  we  have  examined  one  warehouse,  N retailer  pro- 
duction inventory  systems  with  both  continuous  and  dynamic  demand  rates. 
For  the  continuous  demand  rate  case,  we  have  demonstrated  that  a single 
cycle  policy  need  not  be  optimal  for  various  classes  of  production 
policies.  On  the  other  hand,  single  cycle  policies  are  optimal  for 
certain  classes  of  problems  and  the  stationarity  and  renewal  points 
present  in  such  policies  along  with  ease  of  computation  make  them 
especially  appealing  with  regard  to  implementation. 

Further  research  is  indicated  to  determine  bounds  on  the  cost 
difference  between  the  overall  optimal  solution  and  the  best  single  cycle 
policy  for  those  cases  where  it  is  not  the  optimal  solution.  In  addition, 
direct  comparisons  of  the  iterative  procedures  presented  in  this  thesis 
and  the  existing  branch  and  bound  algorithm  (Graves  and  Schwarz,  1977) 
to  determine  single  cycle  policies  should  be  performed.  Particular 
attention  should  be  paid  to  the  relative  effectiveness  of  these  techniques 
for  solving  general  arborescent  problems  through  the  examination  of  many 
one  warehouse,  N retailer  problems. 

Also,  it  was  demonstrated  in  Chapter  2 that  different  one  warehouse, 

N retailer  problems  with  identical  set-up  and  holding  cost  ratios  have 
the  same  single  cycle  solution.  An  interesting  sensitivity  analysis 
question  would  be  how  much  these  ratios  can  be  perturbed  without  changing 
this  identical  solution  property.  Finally,  one  warehouse,  N retailer 
problems  with  less  restrictive  continuous  demand  assumptions  also  warrant 
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investigation.  The  deterministic  single  cycle  solution  may  also  serve  as 
a good  starting  point  in  an  algorithm  to  determine  operating  policies  for 
two-echelon  inventory  systems  with  probabilistic  demands  (Muckstadt,  1977). 

For  the  dynamic  demand  case , we  have  presented  a dynamic  programming 
algorithm  to  determine  optimal  policies  along  with  computational  simpli- 
fications that  can  be  made.  In  addition,  upper  and  lower  bounds  on  the 
optimal  solution  value  were  discussed.  Selected  special  cost  structure 
cases  where  single  facility  techniques  can  be  used  to  determine  optimal 
policies  were  also  introduced. 

Further  research  and  computational  experience  to  determine  the  quality 
of  the  solutions  obtained  by  single  facility  methods  for  general  cost 
structure  one  warehouse,  N retailer  problems  is  indicated.  Similarly, 
the  one  warehouse,  N retailer  problem  with  capacity  constraints  also 
warrants  further  attention.  Dynamic  one  warehouse,  N retailer  problems 
with  stochastic  demands  should  also  be  investigated. 

Finally,  the  techniques  presented  in  this  thesis  should  be  applied 
to  the  other  multi-echelon  structures  both  with  continuous  and  dynamic 
demands.  It  is  possible  that  an  iterative  procedure  similar  to  the  one 
presented  in  this  thesis  can  be  used  to  determine  the  integer  multiples 
of  the  basic  lot  size  in  the  operating  policies  that  have  been  proposed 
for  serial  and  assembly  systems.  In  addition,  the  relation  of  these 
policies  to  overall  optimality  without  the  restrictions  of  stationarity 
and  integer  multiples  can  be  investigated  in  a manner  similar  to  this 
thesis. 

Finally,  the  techniques  developed  in  this  thesis  should  be  applied 
to  determine  operating  policies  for  actual  one  warehouse , N retailer 
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production  inventory  systems.  Such  systems  can  be  found  in  many  industrial 
and  military  situations  and  hopefully  methods  presented  here  can  lead  to 
more  efficient  and  economical  operations. 
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SECURITY  CLASSIFICATION  OF  THIS  PACEflFh«n  D«M  Knfrmd) 


A one  warehouse , N retailer  systen\  is  a special  case  of  the  more 
general  arborescent  production  inventory  system.  We  begin  by  examining 
the  previous  research  on  systems  with  this  arborescent  structure  as 
well  as  reviewing  the  literature  dealing  with  the  serial  and  assembly 
multi-echelon  production  inventory  systems. 

The  continuous  one  warehouse , N retailer  problem  is  then  examined 
in  detail.  The  basic  model  is  introduced  as  well  as  some  previously  solved 
special  cases.  Several  basic  production  policies  that  have  been  suggested 
for  this  system  are  reviewed  and  properties  of  an  overall  optimal  solution 
are  discussed.  We  develop  optimal  and  heuristic  algorithms,  to ^ determine 
the  values  of  the  parameters  in  '‘single  cycle*'  policies  (Schwarz,  1973) 
and  make  comparisons  based  on  both  the  quality  of  solutions  obtained 
and  the  computational  effort.  It  has  been  conjectured  by  Graves  and 
Schwarz  (1977)  that  these  single  cycle  policies  are  optimal  for  certain 
larger  classes  of  production  plans  and  we  demonstrate  that  this  is  not 
the  case.  Finally,  we  distress  the  class  of  multiple  cycling  policies 
and  their  relation  to  optimality/  chseyjssea^j 

We  then  explore  the  dynamic  one  warehouse,  N retailer  problem. 

The  basic  model  is  introduced  arid  both  the  Wagner-Whitin  (1958)  and 
the  Zangwill  ( 1966 )/algorithms /are  reviewed.  Veinott  (1969)  has 
suggested  the  extension  of  Zangwill 's  algorithm  to  the  case  of  N 
retailers,  and  we  present  this  extension /along  with  computational 
simplifications.  Easily  computer  upper  /and  lower  bound/  on  the  optimal 
solution  value  are  developed.  Exact  solution  techniques  for  several 
special  cases  for  the  cost  structure  are  also  discussed.  Finally,  * 

areas  for  future  investigation  are  suggested. 
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